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“La vida en la Tierra ha evolucionado en un mar de campos electromagnéticos naturales. 
Durante el siglo pasado, este ambiente natural ha cambiado en forma significativa con 
la aparición de un vasto y creciente espectro de campos producidos por el hombre. 
Razonando a partir de modelos basados en la termodinámica de equilibrio y en efectos 
térmicos, al principio se pensó que esos campos artificiales eran demasiado débiles 
como para interactuar con los sistemas de moléculas biológicas y por tanto incapaces 
de influir sobre las funciones fisiológicas. Numerosos estudios de laboratorio se 
efectuaron para buscar efectos biológicos a nivel celular y molecular de campos 
electromagnéticos de diferentes frecuencias de oscilación, focalizándose en 
exposiciones que no producen efectos térmicos. Una conclusión que se desprende 
claramente de estos estudios es que muchas de las interacciones observadas no se 
basan en el calentamiento de los tejidos. Los campos electromagnéticos débiles pueden 
modular eventos químicos que ocurren en la superficie de la célula, amplificando los 
efectos de la fijación de hormonas, anticuerpos y neurotransmisores en sus 
correspondientes sitios receptores en las membranas. Los ¡ones calcio juegan un rol 
clave en esta amplificación. Estos estudios fundamentan nuevos conceptos de 
comunicación entre células a través de las barreras constituidas por las membranas 
celulares. Señalan con creciente certidumbre a una organización esencialmente física 
de la materia viva, a un nivel mucho más fino que el correspondiente a la imagen 
estructural y funcional definida en la química de las moléculas.” 


W. Ross Adey, “Biological Effects of Electromagnetic Fields”, Journal of Cellular Biochemistry, 
51: 410-416, 1993. 


“No se ha alcanzado todavía una comprensión a fondo de los efectos biológicos de los 
campos electromagnéticos. El amplio espectro de temas que deben ser estudiados 
incide en el lento progreso hacia ese fin último. Entre las disciplinas involucradas están 
la biología básica, la ciencia médica y la práctica clínica, la ingeniería biológica y la 
ingeniería eléctrica, la química básica y la bioquímica, y la física fundamental y la 
biofísica.” 


Ben Greenbaum, en el capítulo introductorio a la tercera edición del “Handbook of Biological 
Effects of Electromagnetic Fields”, CRC-Taylor 8 Francis, Boca Raton, FL, USA, 2007. 
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(1) Introducción 
En propósito de este informe es presentar algunos aspectos de la conexión entre la 
física del magnetismo y algunos procesos fisiológicos en los organismos. 


En biología se trabaja con un conjunto de niveles estructurales, cada uno de los cuales 
se sustenta en los niveles que lo preceden. Una mirada a los procesos biológicos desde 
la perspectiva del bio-electromagnetismo permite distinguir cuatro niveles: 


En la base, un primer nivel atómico-molecular. Se estudia desde el punto de vista de la 
física cuántica. Generalmente se utiliza la ecuación de Schródinger con términos que 
tienen en cuenta las interacciones eléctrica y magnética en los átomos y entre los 
átomos en las moléculas y entre las moléculas. Es el bien conocido nivel físico 


microscópico. 


Un segundo nivel corresponde a estructuras supramoleculares y sus procesos. 

Se estudian con herramientas de física estadística y de cinética física, teniendo en 
cuenta las interacciones eléctricas y magnéticas, así como los resultados de la 
investigación de las estructuras y los procesos en el nivel previo. Es el nivel físico 
mesoscópico. Se describen conjuntos numerosos a átomos y moléculas, cuyas 
dimensiones son grandes respecto de las distancias entre los átomos, pero son 
pequeñas respecto de las dimensiones características de los patrones de estructura 


(macroscópicos) que aparecen en el siguiente nivel. 


Un tercer nivel, ya macroscópico, corresponde a la célula considerada como una unidad 
de estructura y función. Aquí se incluyen organismos unicelulares, como las bacterias. 
Los procesos a este nivel se estudian, desde un punto de vista fisicoquímico, con la 
ayuda de modelos basados en la termodinámica de los sistemas fuera del equilibrio. 
En las ecuaciones que relacionan los flujos con las fuerzas termodinámicas se incluyen 
los efectos de los campos eléctricos y magnéticos. 

Los parámetros correspondientes los modelos formulados a este nivel se interpretan 
mediante estructuras y procesos pertenecientes a los dos niveles precedentes. 


Finalmente, un cuarto nivel que corresponde al organismo en su totalidad, con sus 
tejidos y órganos. Ahora el estudio se lleva a cabo con las herramientas del 
electromagnetismo clásico de medios continuos. Cuando se construyen modelos 
matemáticos de tejidos desde una perspectiva embriológica o fisiológica, algunas veces 


el nivel celular se trata como si fuera mesoscópico y el cuarto nivel como macroscópico. 
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Las ecuaciones de Maxwell con sus parámetros fenomenológicos (permitividad 
dieléctrica, permeabilidad magnética y conductividad eléctrica) se aplican a modelos de 
medios continuos polifásicos y posiblemente anisótropos. 

En el caso del estudio de la interacción entre un organismo y un campo eléctrico, 
magnético o electromagnético, por lo general no resulta factible tener en cuenta los 
detalles de la distribución espacial de las propiedades eléctricas del tejido. Las regiones 
abarcadas por las corrientes inducidas incluyen varios tipos de tejidos diferentes e 
interfaces entre unos y otros. En estos casos, aún con el apoyo de modelos numéricos, 


es necesario trabajar con valores promediados. 


La aproximación más simple, aunque a menudo insuficiente, es trabajar como si el 
organismo fuera un conductor de volumen homogéneo e isótropo cuya conductividad es 
un valor promedio. No obstante, algunas veces los patrones de variación espacial de la 
conductividad pueden presentar regularidades o tener unas dimensiones tan pequeñas 
que el conductor de volumen de los tejidos puede ser definido macroscópicamente como 
compuesto por partes homogéneas e isótropas. 

El siguiente paso en complejidad es, por tanto, considerar que el conductor de volumen 
de los tejidos se puede dividir en regiones homogéneas e isótropas, pero con 


conductividades diferentes. 


Un modelo más realista se obtiene asumiendo que algunas de estas regiones son 
anisótropas. La Figura 1 corresponde a un corte transversal de un modelo eléctrico de 


tórax humano (Malmivuo y Plonsey, 1995). 


Figura 1: Corte transversal correspondiente a un modelo de tórax polifásico, heterogéneo 
(diferentes resistividades en diferentes regiones) y anisótropo (diferentes resistividades 
según la dirección respecto a las fibras musculares esqueléticas o miocárdicas). 
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La resistividad disminuye a medida que los colores pasan de amarillo claro (que 
corresponde a huesos y cartílagos) al rojo (que corresponde a la sangre en los 
ventrículos cardíacos y en la aorta). 

Algunas veces se pueden introducir simplificaciones adicionales, dependiendo del 
problema considerado. 

Si lo que interesa son los campos en las cavidades del corazón, se puede sustituir el 
modelo más completo de la figura precedente por un primer modelo más simple en el 
cual, además de las cavidades cardíacas, se tienen en cuenta dos fases homogéneas, 
una correspondiente a los pulmones y la otra al resto de los tejidos (Figura 2 A, Malmivuo 


y Plonsey, 1995). 


Figura 2 A: Simplificación del modelo de la Figura 1. 


En esta figura las resistividades de cada fase aparecen en Qxm. Cada una de las tres 


fases se caracteriza por la región que ocupa y por su conductividad. 


Un segundo modelo, aún más simplificado, se puede ver en corte en la Figura 2 B 
(Malmivuo y Plonsey, 1995). Este modelo distingue solamente dos fases: una fase que 


representa a las cavidades cardíacas y otra fase que representa al resto de los tejidos. 


> 


Figura 2 B. Simplificación del modelo de la figura 2 A. 


4 


Septiembre de 2020 


Considerados globalmente, los tejidos biológicos y e incluso sus células, se comportan 
como medios diamagnéticos. Los eritrocitos desoxigenados son una excepción, pues 
se comportan como paramagnéticos. Además, en bacterias y tejidos biológicos de los 
animales y las plantas puede haber partículas ferrimagnéticas (formadas de magnetita). 


La exposición a campos ambientales, a campos asociados a instalaciones industriales, 
a campos muy intensos en el diagnóstico médico (imágenes de resonancia magnética 
nuclear) y en el tratamiento médico (estimulación magnética funcional) es cada vez más 
frecuente. Esto ha aumentado el interés por los mecanismos de interacción de los 
campos con los tejidos biológicos y la preocupación por establecer umbrales de 
exposición que no deben ser superados. 


En este trabajo el énfasis se pone en estudiar de algunos de los fenómenos que 
involucran el campo magnético en las estructuras y procesos correspondientes a los 
últimos dos niveles de la jerarquía. En ocasiones haremos referencia a resultados 


obtenidos en la investigación de los dos primeros niveles. 


En particular consideraremos algunos de los fundamentos físicos de la interacción de 
los campos magnéticos ambientales con los organismos: 

-Orientación de partículas en campos magnéticos estáticos. 

-Traslación de partículas en campos magnéticos no homogéneos. 

-Aparición de fuerzas sobre partículas cargadas en movimiento en campos magnéticos. 
-Inducción de corrientes eléctricas en el conductor de volumen de los tejidos debidas a 
variaciones del campo magnético en el tiempo o al movimiento de los organismos 


inmersos en campos magnéticos estáticos que presentan variaciones espaciales. 


Para el caso de campos estáticos homogéneos, estudiamos brevemente la orientación 
y deformación de células y de estructuras supramoleculares debido a fuerzas y a torques 
magnéticos. En tejidos isótropos esas fuerzas de volumen (fuerzas de traslación) 
aparecen cuando la susceptibilidad magnética del tejido no es uniforme. 

Los torques pueden actuar sobre las partículas de magnetita mencionadas previamente. 
También aparecen torques en campos homogéneos debidos a anisotropías en la 
susceptibilidad magnética de estructuras supramoleculares diamagnéticas, como las 
membranas celulares y los cromosomas. 

Para campos magnéticos estáticos no homogéneos, consideramos la aparición de 
fuerzas magnéticas y sus efectos. 

Para campos dinámicos, revisamos la generación de corrientes eléctricas inducidas en 


el conductor de volumen de los tejidos, que cuando disipan energía siempre producen 
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calor. Si poseen las características adecuadas esas corrientes inducidas pueden 


estimular a nervios, a músculos esqueléticos e inclusive al músculo cardíaco. 


A diferencia del campo eléctrico, el campo magnético no es atenuado en la interfaz entre 
el aire (o el agua) y el organismo, debido a que las diferencias en las permeabilidades 
magnéticas de estos medios difieren muy poco entre sí, al igual que las diferencias en 
la permeabilidad que se encuentran entre diferentes tejidos. 

Esto es así exceptuando las partículas ferri-magnéticas de dimensiones del orden de 
decenas de nanómetros o de micrómetros, que pueden aparecer en el interior de 
algunas bacterias y en tejidos de animales. Estas partículas pueden tener un rol 
fisiológico, o pueden ser el resultado de la contaminación de tejidos (pulmones, hígado) 
por inhalación o ingesta. * 

Las mediciones de la susceptibilidad magnética del hígado y de los pulmones 
suministran indicadores del grado de esa contaminación. 

Ubicando el tórax de una persona, cuyos pulmones se han contaminado por inhalación 
de polvo con contenido ferrimagnético, en un campo magnético constante lo bastante 
intenso durante algunos segundos, las partículas de magnetita se orientan 
paralelamente al campo externo. Cuando éste desaparece, el campo residual debido a 
los momentos magnéticos sufre un proceso de relajación. La evolución del campo 
residual se puede medir y se obtiene lo que se denomina magneto-neumografía. 

Este procedimiento suministra información sobre la contaminación en los alvéolos, 
sobre parámetros viscoelásticos celulares y sobre el movimiento de células del sistema 
inmunitario que capturan y eliminan material extraño al tejido. 

Aunque las diferencias en la susceptibilidad magnética entre los diferentes tejidos 
diamagnéticos son mucho menores, la diferencia de susceptibilidad magnética entre el 
músculo cardíaco y la sangre en las cavidades ventriculares es suficiente para permitir 


la determinación de volúmenes. 


El campo magnético puede actuar sobre momentos magnéticos a nivel nuclear 
(interacción empleada en los métodos basados en la resonancia magnética nuclear) y 
atómico (efecto Zeeman extrínseco). Además, el campo puede interactuar a niveles 
moleculares y supramoleculares, y sobre las corrientes eléctricas circulantes en los 
tejidos biológicos, relacionadas con algunos los procesos fisiológicos. 

Si los procesos en un organismo se ven influidos por un campo magnético, ese 


organismo puede ser considerado como un detector de ese campo. Si así se lo 


1 Las partículas que poseen un rol fisiológico difieren desde el punto de vista 
cristalográfico de las partículas provenientes de la contaminación. 
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considera, un campo magnético externo podría pensarse como una señal. Para ser 
detectada, debe superar el nivel de ruido. Este es otro de los temas que abordaremos 
en lo que sigue. 


También estudiamos los fundamentos físicos de la generación de campos magnéticos 
por esos mismos organismos y los fundamentos de la medición de esos campos 
mediante instrumentos externos, como el SQUID. 

En medicina el magneto-cardiograma, el magneto-encefalograma y la magneto- 
miografía suministran información complementaria a la obtenida, respectivamente, del 


electrocardiograma, el electroencefalograma y la electromiografía. 


En general no es posible desvincular totalmente los campos eléctricos que aparecen en 
el conductor de volumen electrolítico formado por los tejidos biológicos, de los campos 
magnéticos ambientales o producidos por los organismos. 

Se pueden distinguir dos fuentes principales de los campos magnéticos producidos por 
los organismos. 

Una de estas fuentes son las corrientes eléctricas que acompañan a los procesos 
fisiológicos en el corazón, el sistema nervioso y los músculos, o menos frecuentemente 
y con menor importancia, las corrientes eléctricas asociadas al campo eléctrico 
generado cuando el organismo se mueve respecto de un campo magnético externo. 
La otra fuente del campo magnético producido en los organismos son los cuerpos 
magnetizados presentes en los tejidos biológicos o en el interior de las bacterias. 


Cuando se tienen en cuenta las interacciones de los tejidos biológicos con campos 
externos variables con el tiempo, en condiciones casi-estacionarias (es decir, cuando se 
puede despreciar la corriente de desplazamiento en el interior de los tejidos) un campo 
magnético variable local induce un campo eléctrico local variable. 

Pero el campo eléctrico variable por sí mismo (es decir, sin la intermediación de las 
corrientes eléctricas asociadas a esos campos) no genera un campo magnético en estas 
condiciones. 

Pero a frecuencias lo bastante elevadas como para no poder despreciar la corriente de 
desplazamiento, el campo eléctrico local variable induce directamente un campo 


magnético local y se transforma entonces en una fuente de campo magnético. 


Puesto que los seres humanos hemos estado expuestos durante eones a los campos 
eléctricos y magnéticos naturales, presenta interés revisar algunos aspectos 
relacionados con las amplitudes y frecuencias de estos campos para poder compararlos 
con los campos originados en las actividades propias de nuestra Era Tecnológica 
(Suárez-Antola, 2020). Comenzaremos entonces por este tema. 
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(2) Campos geomagnéticos y geo-eléctricos estáticos. Orientación de 
bacterias en el campo magnético terrestre. 


En la superficie de nuestro planeta, las plantas, los animales y los seres humanos se 
encuentran con campos estáticos de inducción magnética B de magnitudes 


comprendidas entre 33 y 67 1T y campos eléctricos también estáticos E con magnitudes 
en un entorno de los 100 V/m, pudiendo alcanzar los 100 kV/m bajo una nube de 
tormenta. Estos campos estáticos no parecen estar implicados en efectos biológicos. 
El campo geo-eléctrico estático no parece producirlos en ausencia de vías de descarga 
a tierra que permitan la instalación de corrientes eléctricas significativas en el conductor 
de volumen de los tejidos. 

Pese a que en la mayor parte de los organismos el campo geomagnético estático no 
parece producir efectos biológicos en ausencia de campos oscilantes que se le 
superpongan, en algunos organismos el campo geomagnético estático se relaciona con 
la posibilidad de un cierto grado de orientación espacial. 

Algunas bacterias que viven en medios acuático (lagos o mares) y son anaerobias, 
poseen en su interior un dipolo magnético m en forma de barra de varios micrómetros 
de longitud formada por pequeños cuerpos (magnetosomas) de decenas o centenares 
de nanómetros, magnetizados y adheridos entre sí, con la misma orientación de su 
momento magnético. 

Los magnetosomas están formados por una partícula de magnetita (ferrimagnética) ? 
rodeada por una membrana. 

Cada partícula constituye por sí misma un único dominio magnético. 

El campo magnético terrestre produce un torque 7,,, = mi x By sobre ese dipolo. 

Este torque gira la bacteria tendiendo a alinearla con el campo magnético terrestre. 

La bacteria alineada con el campo se desplaza activamente utilizando sus cilias o su 
flagelo. * 

Si la magnetización de la barra de magnetosomas es la apropiada, el microorganismo 
se aleja de la superficie rica en oxígeno (donde no puede vivir) y se dirige a las regiones 
profundas, pobres en oxígeno, donde puede vivir y reproducirse. 

Cuando una de estas bacterias se reproduce, la mitad de la barra de magnetosomas 
pasa cada célula bacteriana hija. 


2 En biología y medicina a menudo se hace referencia a ellos como si fueran 
ferromagnéticos. 


3 El movimiento en sí resulta de la actividad de la bacteria, pero la orientación en el campo 
magnético externo es puramente pasiva. 
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Estas sintetizan nuevos magnetosomas que se añaden a los recibidos de la célula 
bacteriana madre hasta que las longitudes de las barras alcancen valores que aseguren 
torques adecuados para girar la célula bacteriana en el campo magnético terrestre. 

La Figura 3 (Glaser, 2012) muestra un esbozo de algunas líneas del campo magnético 
de la Tierra. 

A la derecha se ve la inclinación del campo magnético en dos localidades. 

En la figura se puede ver la descomposición del campo en dos: una componente vertical 
(que es máxima en las cercanías de los polos y mínima en las proximidades del ecuador) 
y una componente horizontal (que es máxima en las proximidades del ecuador y mínima 


en las cercanías de los polos). 


Figura 3 Esquema del campo magnético de la Tierra con un detalle de la descomposición 
del campo en una componente vertical y una componente horizontal en dos localidades 


señaladas con círculos en la superficie terrestre. 


Desde hace mucho tiempo se dispone de evidencia acerca de la orientación magnética 
de animales (insectos, aves, peces y mamíferos) mediante diversos tipos de receptores 
geomagnéticos, pero todavía no hay suficiente consenso entre los investigadores 


acerca de los mecanismos involucrados en cada caso. 


(3) Campos geomagnéticos y geo-eléctricos variables. Comparación con los 
campos generados por las líneas de transmisión de potencia. 

Como los campos geo-eléctricos y geomagnéticos presentan fluctuaciones asociadas a 

las tormentas atmosféricas y a las perturbaciones en las cargas de la ¡onósfera, todos 
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los organismos biológicos se encuentran inmersos en campos variables que se 
superponen a los campos estáticos. 

En el caso del campo eléctrico, aparecen fluctuaciones naturales de 0.01 V/m a 
frecuencias hasta de 10 Hz. 

Esas fluctuaciones decrecen hasta alcanzar valores cien veces menores, comprendidos 
entre 0.1 y 0.6 mV/m a frecuencias de 50 Hz. 

A frecuencias mayores a algunos centenares de Hz la amplitud de oscilación del campo 


eléctrico, superpuesta a su valor estático, es aún menor. 


Resulta interesante comparar estas amplitudes naturales de oscilación del campo 
eléctrico a 50 Hz con las amplitudes del campo eléctrico oscilante a esa misma 
frecuencia debajo de una línea de transmisión de potencia trifásica, que puede alcanzar 
unas decenas de kV/m y son capaces de encender un tubo fluorescente que 
resplandece sin alambres de conexión cuando se lo mantiene en el aire debajo de la 


línea. La relación de amplitudes es de 1 a 10”. 


En el caso del campo geomagnético, las amplitudes de las oscilaciones a 1 Hz son de 
1 nT, mientras que a frecuencias comprendidas entre 10 y 100 Hz dichas amplitudes 
son de 10 pT. 

Si bien durante una tormenta magnética la intensidad puede aumentar en 0.5 u'T, este 
aumento se produce en minutos u horas. 

Debido a las influencias del sol y de la luna sobre las corrientes ¡ónicas en la alta 
atmósfera, aparecen variaciones diurnas, con amplitudes de aproximadamente 30 nT. 
Además, superpuestos a esos campos lentamente variables asociados a eventos 
atmosféricos irregulares, aparecen otros campos, débiles, generados por un fenómeno 
de resonancia en la cavidad comprendida entre la superficie terrestre y la capa más baja 
de la ionósfera. La energía de estos campos proviene de la descarga de rayos. 

La formación de campos periódicos con frecuencias de aproximadamente 8 Hz, 14.1 
Hz, 20.3 Hz y 32.5 Hz se debe a los fenómenos de propagación y reflexión en esa 
cavidad resonante. Este fenómeno se conoce como resonancia de Schumann. 

Las amplitudes de las distintas componentes disminuyen al aumentar la frecuencia, pero 
en general aún las más bajas no superan los 20 nT (generalmente son bastante 
menores). 

Es un hecho interesante la variación en las frecuencias de Schumann provocadas por 
modificaciones en el espesor de la cavidad resonante. Ese espesor se modifica debido 
los efectos gravitacionales del sol, la luna y de algunos planetas sobre la ¡onósfera. 
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Resulta asimismo interesante comparar la amplitud de oscilación de 10 pT del campo 
magnético a 50 Hz con los más de 0.1 uT de amplitud de oscilación a esa misma 
frecuencia, que se pueden medir debajo de algunas líneas de transmisión de potencia. 
La relación de amplitudes es de 1 a 10*. 


En relación con estas comparaciones, es conveniente llamar la atención sobre el 
apantallamiento que sufren los campos eléctricos ambientales, naturales y artificiales, 
por árboles, edificios y otros objetos. 

Ese apantallamiento en general no se manifiesta en el caso del campo magnético 
ambiental, natural o artificial, que disminuye fundamentalmente cuando aumenta la 
distancia entre su fuente y el punto de observación. 

Por otra parte, el campo eléctrico en un cierto punto situado debajo de una línea de 
transmisión de potencia depende del voltaje con el cual se transporta la energía, cuya 
amplitud es constante. En cambio, el campo magnético en ese mismo punto varía 


cuando la corriente que circula por la línea se modifica durante el curso de un día. 


Como los campos eléctricos y magnéticos artificiales son mucho mayores que los 
campos naturales que oscilan a esas mismas frecuencias y que nos han acompañado 
durante nuestra evolución biológica como especie, parece que debiera ser de interés el 
estudio de sus posibles efectos sobre nuestro organismo. 


(4) Ruidos, campos y contaminaciones. 
La mayor parte de las personas que investigan los efectos biológicos de los campos 
electromagnéticos asume que los campos y corrientes asociados con el funcionamiento 
normal del organismo no deben producir ellos mismos (ni deben actuar como 
promotores de) daños significativos a nivel de las células, ni sobre sus orgánulos 
funcionales ni sobre el ADN. 
Entonces parece que no se deberían esperar efectos perjudiciales provocados por las 
fluctuaciones aleatorias naturales de los campos eléctricos y magnéticos a nivel de las 
membranas celulares. 
Lo mismo cabe esperar en el caso de los campos magnéticos generados por las 
corrientes eléctricas que normalmente circulan por el organismo, asociadas a las 
actividades del sistema nervioso central, del sistema neuromuscular y del corazón. 
Como dijimos en la introducción, si un organismo se ve afectado de cualquier manera 
por un campo externo, parece que se lo podría considerar como un detector natural para 
ese campo. 
El campo externo sería una señal que recibe el detector, señal que para ser detectada 
debe superar al ruido. 
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Suponiendo que esto sea así, la intensidad de estos campos de ruido suministra una 
escala natural para comparar las intensidades de los campos intracelulares o en el 
interior de las membranas, generados por la exposición del organismo a campos 
externos. Si las intensidades de estos últimos campos son mucho menores o a lo sumo 
del mismo orden de magnitud que los campos de ruido, parecería que no deberían 
producir efectos significativos. 

Los magnetosomas pueden resultar de importancia en relación con posibles efectos 
sobre las estructuras supramoleculares de las células, debidos a un campo magnético 
externo. 

En 1992, Kirschvink y otros (1992) informaron que el cerebro humano contiene varios 
millones de micropartículas ferrimagnéticas por gramo de tejido: los magnetosomas que 
ya mencionamos a propósito de la orientación de algunas bacterias en el campo 
magnético de la Tierra. Posteriormente Kobayashi y otros (1995) descubrieron que la 
contaminación con magnetosomas puede afectar los experimentos bioeléctricos y 
biomagnéticos realizados con cultivos de tejidos cuyas células normalmente no 
contienen magnetosomas. 

Esta contaminación puede deberse a los materiales plásticos previamente esterilizados 
que se utilizan en experimentos de laboratorio con cultivos de tejidos: contienen 
partículas ferrimagnéticas de menos de 100 nm, que son incorporadas fácilmente por 
las células de la serie blanca de la sangre. 


Para estimar el cociente señal-ruido para uno de estos pequeños magnetosomas en un 
campo magnético estático, se puede asumir una situación de equilibrio local, como se 
hace en la termodinámica de procesos irreversibles. 

Esta suposición puede hacerse para la mayor parte de las redes de reacciones 
bioquímicas y procesos de transporte en el interior celular y en las membranas 
biológicas, aunque no en todos los casos. En particular, puede que no sea posible hacer 
esta suposición durante algunas de las bruscas transiciones de fase fuera del equilibrio 
que se producen a nivel celular. 

Admitiendo equilibrio local y la posibilidad de una descripción aplicando la física 
estadística clásica, la densidad de probabilidad f;, (9) correspondiente al ángulo 9 que 


el momento magnético m de un magnetosoma forma con el campo magnético terrestre 


Br viene dada por: 


ON (27 sen 0 exp E) / (4 (22) senh A ] (4.1) 


T 
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Esta expresión es formalmente la misma que se aplicó originalmente en un análisis 
microscópico del paramagnetismo basado en la física estadística clásica (Sommerfeld, 
1956). Ahora la estamos aplicando en un análisis a nivel mesoscópico. 


m Brycosó 


ha T | se reduce a la 


En ausencia de campo externo el factor de Boltzmann exp | 
unidad para todos los ángulos posibles y la densidad de probabilidad de que el momento 


is ; 1 1 da 
magnético forme un ángulo 8 con el campo se reduce a pa 21 sen0 = Po donde 


o = 211 senO9 Aquí dQ = 271 sen O d6 es un elemento de ángulo sólido y 47 es el 


ángulo sólido total, de modo que - es la probabilidad de que, en ausencia de campo 


externo, el ángulo entre el momento magnético y el campo se encuentre entre O y 0 + 
d0. En la Figura 4 siguiente se puede ver un esquema un elemento de ángulo sólido 
dQ = 21 sen 6 d0 correspondiente al caso en el que las velocidades de desplazamiento 


de la bacteria forman un ángulo 6 con la dirección y sentido del campo magnético. 


a 2 
a me > Ti - La a 
ALIS 
O AL "dol 


ns 
a 


Figura 4. El eje vertical corresponde a la dirección del campo magnético externo. 
La banda rayada representa el elemento de ángulo sólido dí. = 27. sen O d9 
Sobre esa banda se muestra un vector de momento magnético. 


La componente del momento magnético del magnetosoma paralela al campo es 
mcosóO . 


El promedio de esta componente paralela al campo viene dado por la fórmula: * 


4 Utilizamos una barra por encima de la magnitud que se promedia para indicar un 
promedio tomado en el sentido de la física estadística. El símbolo (..., ...) se reserva para 
ser utilizado en desarrollos que involucran mecánica cuántica. 
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(4.2) 


m, = mcos0 = So fn(0) (mcos 0) de = m [cotgh E =- Ae] 


kg T m Br 
La función £L(x) = cotgh(x) =: fue introducida por Langevin. * Si x es pequeño 
respecto de la unidad o a lo sumo x 3 1, L£(x) se comporta como =x mientras que si 
es grande respecto de la unidad L(x) se encuentra muy próxima a 1, comportándose 
como 1 — - : 


Así, cuando x = 10 la función de Langevin toma un valor próximo a 0.9, de modo que 
la dirección del momento magnético se encuentra en promedio muy cercana a la 
dirección del campo externo. 

El campo magnético tiende a alinear mientras que las interacciones térmicas con las 


partículas del medio tienden a desalinear el momento magnético del magnetosoma. 


ñ E m 4 B Y 
En consecuencia, el número no-dimensionado x = A se puede interpretar como una 
B 


medida del cociente señal magnética/ruido térmico, donde m B es la energía máxima 
de interacción magnética entre la partícula magnetizada y el campo magnético externo 
y kg T es una estimación de la energía del ruido térmico local. 

Cuando ese número posee un orden de magnitud de la unidad o menor, el ruido térmico 
tiende a romper la alineación del momento magnético con el campo externo y 


prepondera el desorden molecular sobre el orden asociado a la alineación. 


; mB A s e 
Por el contrario, cuando ES de un orden superior a la unidad, la señal supera al 
B 


ruido y la alineación se puede estabilizar, aunque con fluctuaciones tanto más 
frecuentes cuanto menor sea el orden numérico del número no-dimensionado. 
Supongamos que un magnetosoma se encuentra en un campo magnético terrestre de 
5x10* T, actuando a una temperatura local de 300 *K. 

En general las partículas magnéticas que no se deben a la contaminación y parecen 
cumplir un rol fisiológico poseen un diámetro próximo a los 50 nm, alcanzando algunas 
veces los 100 nm. Si su diámetro es demasiado grande, se podrían formar varios 
dominios con diferentes orientaciones magnéticas. En ese caso el momento magnético 
se vería disminuido respecto del caso en el cual un único dominio abarca la totalidad de 
la partícula. Si su diámetro es demasiado pequeño, las fluctuaciones térmicas en el 
medio podrían destruir el alineamiento del momento magnético de la partícula con el 


campo externo. 


5 Langevin obtuvo la función L(x) en una investigación, empleando herramientas de física 
estadística clásica, de las propiedades dieléctricas de un gas ideal, cuyas moléculas 
poseen momentos dipolares permanentes, en un campo externo. Posteriormente extendió 
este enfoque al paramagnetismo (Sommerfeld, 1956). 
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Si el diámetro de un magnetosoma es de 50 nm cabe esperar un momento magnético 


de 6.4x10"'” Axm? y en este caso a adopta como valor 0.77. 
B 


Si su diámetro es de 100 nm y está magnetizado formando un dominio único, su 


momento magnético es de 2x 10? Axm? y ahora =- vale 24. 
B 


Para magnetosomas con diámetros crecientes entre los 50 nm y los 100 nm, el número 
mB ] A 
e y toma todos los valores intermedios entre 0.77 y 24. 
B 
Así pues, debido a que puede dar lugar a una elevada relación señal-ruido en el campo 
magnético terrestre o en campos externos más intensos, la existencia de 
magnetosomas, integrados en estructuras supramoleculares y eventualmente actuando 
en forma sinérgica, podría tener como consecuencia la aparición de algún efecto 


biológico, como la apertura de canales iónicos en las membranas celulares. 


(5) Un modelo simplificado de efecto fisiológico del campo magnético. 
Consideremos entonces un modelo de apertura de poros en una membrana biológica, 
asociada a la rotación de un magnetosoma esférico de radio * a y de momento 
magnético m = 2 x 10% a? Axm?, en un campo magnético variable que se superpone 
al campo magnético terrestre. 

Supongamos que el magnetosoma se halla integrado en la membrana, adherido a la 
compuerta de un canal iónico. 

Además, supongamos que cuando el magnetosoma se encuentra en su posición de 
equilibrio alineado con el campo magnético terrestre, una compuerta, también en su 
posición de equilibrio, obtura el canal. 

Se añade un campo externo que varía con el tiempo B,(t) perpendicular al campo 
magnético terrestre By (que se considera constante) que tiende a rotar al magnetosoma 
en el plano determinado por ambos campos, un cierto ángulo 6(t) respecto de su 
posición de equilibrio alineado con B+ . 

Cuando el magnetosoma rota arrastra tras de sí a la compuerta que rota a su vez un 
mismo ángulo 6(t) y como consecuencia el canal se abre. 

Cuando la apertura alcanza un umbral 6,, , se desencadena un proceso regenerativo 


que permite que un pulso de corriente de ¡ones de cierto tipo atraviese la membrana. 


$ El radio a se debe expresar en m. 
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A medida que la compuerta se aleja de su posición de equilibrio, aparece un torque 
T¿( 0) que se opone a la apertura. 

Para los fines de construir un modelo matemático simple asumiremos que lo podemos 
linealizar introduciendo una constante elástica k. : Te = —k¿0 

A este torque se le suma un torque restaurador debido al campo magnético terrestre 


Tr =—m By seng , un torque perturbador 7, =mB,(t)cosg9 debido al campo 

externo perturbador, un torque Tpr(t) relacionado con el ruido térmico y un torque 
d : 

restaurador T, = —QA Era) debido a fuerzas viscosas. 


Si la viscosidad del medio es u resulta: a = 8r y a? 

Los valores medidos para 87 y se ubican generalmente entre 0.003 y 0.015 Nxs/m? lo 
que sugiere tomar como valor representativo, para =0,0035 Nxs/m?: «a = 0.009 a* 
El momento de inercia / respecto del eje de rotación que pasa por el centro de la esfera, 


suponiendo que el material ferrimagnético posee una densidad uniforme ” p viene dado 
8 5d mia iZ 
por: /= TP a? Entonces la ecuación del movimiento de rotación del magnetosoma 


adherido a la compuerta de un canal que atraviesa la membrana celular se puede 
escribir así (Kirschvink y otros, 1992): 


d? 
dt? 


I 0=-—a 8 —k¿0—m By seno +m B,(t) cos 0 + Tp[t) (5.1) 


Para concretar, consideremos que el campo externo oscila con una frecuencia angular 
WwW y una amplitud Bo: Bp(t) = By cos w t (5.2) 
Para simplificar el modelo, en una primera aproximación no tomaremos en cuenta el 


torque de ruido Tr(t) , aproximaremos sen0 = 6 y cosO = 1. 


, y E j WI wa? , 
Si el número sin dimensiones * a (1 < 1 , el torque viscoso es lo bastante 
intenso como para no tener en cuenta en la ecuación el término debido a la aceleración 
angular. 
Como para una partícula de magnetita de 50 nm de densidad 5240 kg/m? inmersa en 


2 
un medio con una viscosidad =0,0035 Nxs/m?, (2) => es del orden de 107*% x w, y 


como las frecuencias angulares de los campos oscilantes considerados son de un orden 


7 La magnetita (Fez04) posee una densidad de 5240 kg/m?. 


$ Este número se obtiene sustituyendo una solución armónica de frecuencia angular w en 
E d a? pr es d e z 
el término inercial / mE 0 y en el término disipativo a PF 0 , dividiendo luego la amplitud 


del resultado obtenido para el término inercial por la amplitud del resultado obtenido para 
el término disipativo. 
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numérico inferior a 10% Hz, en lo que sigue asumiremos el dominio de la disipación sobre 
la inercia. 


En ese caso la ecuación del movimiento se reduce a: 


cos(ow t) (5.3) 


(k¿+m Br) m Bo 


a 


O = 


d 

E 0+ 

La solución estacionaria de esta última ecuación es: 0(t) = O, cos(w t — p(w)) con 
am w m Bo Bo 


SONS a A y 0. 


k¿+m Br =) Br 3 Ja? w2+(k¿+m Br)? ley ER En) 


. a ' A ñ 
Como tanto a« como m son proporcionales a a?, E resulta independiente del radio del 


magnetosoma. 

Cuando 0,» > 0,, la compuerta rota, alcanza y supera el umbral. El pulso iónico a través 
del canal se produce: el campo magnético externo habrá dado lugar a un efecto 
biológico. 


(6) Influencia de las fluctuaciones estadísticas sobre la orientación de 
partículas magnetizadas. 


Supongamos que el campo B,,(t) desaparece y una partícula magnetizada permanece 
en el campo magnético de la tierra. 

En equilibrio local se puede aplicar la distribución de Boltzmann para una energía de 
interacción U(O) = —mBycosó con0<0<T. 


Para ángulos lo bastante pequeños, si aproximamos cos O por la fórmula aproximada 


mBr 


pg? 


1 =7 0? obtenemos la aproximación U(O) = —m Br + 


Si admitimos trabajar con esta última expresión de la energía potencial sin limitaciones, 
se puede aplicar el teorema de equipartición de la energía de la física estadística para 


obtener el valor medio cuadrático 0,,. del ángulo entre el momento magnético y el 


lez kg T 
campo en un medio de temperatura local T : Omc = 10? = EE 
T 


Asumiendo un momento magnético m = 2 x 10% a? Axm? y un campo terrestre de 


magnitud By =5x 107” T resulta: Oo vs E (6.1) 


Esta expresión del valor medio cuadrático del ángulo entre m y By disminuye como el 


recíproco del cubo del diámetro de la partícula. 
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La razón por la cual 06,». no está acotado superiormente cuando el diámetro tiende a 
cero es la validez de la aproximación cosO = 1 — > 9? que permitió aplicar el teorema 


de equipartición de la energía. 


Solo si O es lo bastante pequeño la aproximación resulta aceptable. 


Para un magnetosoma ligado a la compuerta de un canal iónico, la energía potencial en 


el campo magnético terrestre se puede aproximar así, añadiendo la energía potencial 


m 


Er g2 4 lc g2 
2 


elástica de la compuerta rotada: U(0) = —m By + A 


AZ kg T 
Entonces: Omc = E = ni (6.2) 


Asumiendo nuevamente un momento magnético m = 2 x 10% a Axm? y un campo 


terrestre de magnitud By =5 Xx 107? T resulta: Omc = 10007 mn 


iz TA ; dd 1 Ae ñ 
También para este se utilizó la aproximación cos0 =1 — a 9? con las limitaciones que 


conlleva. 


Cuando es superado el umbral de rotación 0,, , se desencadena la apertura del canal. 
Si f. (0) es la densidad de probabilidad de que el momento magnético, ligado a la 
compuerta, en presencia del campo geomagnético, rote un ángulo € debido a una 
fluctuación aleatoria, la probabilidad de apertura espontánea del canal viene dada por: 


Se fin (0) d0 . Si se efectúa la aproximación cos 9 = 1 =5 0? se puede utilizar, en 


lugar de f,, (9) , una distribución de Boltzmann, basada en el factor exponencial 


La integral se toma en ese caso entre 0,, y +0. 


(7) Aplicaciones de las fuerzas de traslación debidas a campos magnéticos no 
homogéneos a la separación de células. 
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Cuando una partícula magnetizada se encuentre en un campo de inducción magnética 


que, para un mismo instante de tiempo, varía de un punto a otro del espacio, sufre una 


> 


fuerza de traslación: ? En = (m.V) B (7.1) 


Esta expresión se reduce a EF, = vV(m.B ) cuando el momento magnético es 


intrínseco. 
Ñ Ñ Eo ñ > 1 
Cuando el medio es isotrópico y lineal: F,, = V E m.B ) 
Como es usual, el operador V(  ) representa el gradiente, el operador Vx representa 
el rotor del campo al que se aplica, y el operador V- representa su divergencia. 


En coordenadas cartesianas: m.V( )=my =( )+ m, Al )+m, =( ) 


La expresión de la fuerza en un campo magnético no homogéneo se puede deducir a 
partir de un lazo de corriente equivalente al momento magnético. 


Cuando el campo es meridiano-plano, en coordenadas cilíndricas se tiene: 


B(r,z) = B,(r,z) 8, + B, (r,z) €, (7.2) 
Si el momento magnético se encuentra ubicado sobre el eje z de simetría del campo: 
m =m e, (7.3) 


Teniendo en cuenta que sobre el eje de simetría r = 0 la expresión de la fuerza 


magnética se reduce a: En (0, z) = FEn,z(0,Z) €, = (m —B,(0, 2)) €, (7.4) 


2 La expresión F,, = (m.V) B para la fuerza se puede deducir a partir de la fórmula de 
la fuerza que actúa sobre un lazo delgado a través del cual circula una corriente ¿ cuando 


se encuentra en un campo magnético externo B: F,, = $ dí xB A partir del 


teorema de Stokes se puede demostrar la igualdad: $dr x B = ff(dS x V) x B 
donde la integral de superficie se toma sobre una superficie abierta, orientada y lo 
bastante regular, que tiene como borde el lazo de corriente. En los elementos diferenciales 


vectoriales dí = (ds) É (longitud) y dS = (dS) Tí (área), como es usual, í es el versor 
tangente en cada punto del lazo y ni es el versor normal en cada punto de la superficie 


bordeada por el lazo. Como en todas partes se verifica V : B=0 y como en puntos a la 
vez exteriores al lazo de corriente y a las tuentes del campo se verifica V Xx B =0 resulta 
que (as Xx v) x B = (as ] v)B . Como consecuencia de esta última igualdad: F= 
i$bdr xB=if (as : v)B Tomando el lazo lo bastante pequeño es posible sustituir 
¿ff(dS-V)B por (( ff ds) - v) B de modo que: F,=i$dr xB= 
(i([f ds) - V) B_ Pero i ff dS =mi es el momento magnético del lazo, entonces 


resulta finalmente: F,, =(m.V)B Puesto que, desde la perspectiva de la 
electrodinámica de medios continuos, todo momento magnético equivale a un lazo de 
corriente adecuadamente definido, se obtiene el resultado buscado (Kompaneyets, 1978, 
volumen 2; Landau y Lifshitz, 1960, volumen 8). 
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Una aplicación interesante de las fuerzas de traslación en campos no homogéneos es 
la separación magnética de células, utilizando partículas super-paramagnéticas 
(partículas paramagnéticas cuya susceptibilidad magnética relativa *” alcanzar valores 
muy grandes respecto de 1 y se comportan en forma lineal en un campo externo). 
Esas partículas exhiben sus propiedades magnéticas solo cuando se instala un campo 
magnético. 

Se fijan partículas super-paramagnéticas de unos 50 nm de diámetro en puntos 
adecuados de un anticuerpo. 

Este anticuerpo es capaz de unirse en forma específica al tipo de célula que interesa 
separar. 

Con este procedimiento se pueden adherir un hasta centenar de partículas a cada 
célula. 

Se pone la muestra en un campo magnético no homogéneo del orden de 1 T con un 
gradiente del orden de 10 T/m. 

Las células que interesan poseen, adheridas a su superficie, combinaciones de 
anticuerpos con partículas super-paramagnéticas, debido a lo cual se pueden separar 
de las otras células. 


Otra aplicación de los campos no homogéneos es la entrega de fármacos ligados a 
micro-transportadores magnéticos en ciertos sitios del organismo. 

La agitación magnética se emplea para modular la liberación de macromoléculas, 
contenidas en el espacio de poros de partículas de polímero que poseen inclusiones 


magnetizadas. 


(8) Torques originados en la anisotropía en las propiedades diamagnéticas de 
las estructuras moleculares y supramoleculares. Fuerzas de traslación 
asociadas a variaciones espaciales de la susceptibilidad magnética. 


Cuando los materiales biológicos diamagnéticos y anisótropos, como la fibrina y el 
colágeno, o las membranas y los cromosomas, se exponen a campos magnéticos 
homogéneos, sufren un torque que tiende a rotarlos hacia la dirección determinada por 
sus propiedades de anisotropía. 


10 Cuando la relación entre el campo de inducción B' y el campo magnético H_ adopta la 


forma B = Hp (H + M) la relación M = y H para un medio isótropo involucra la 
susceptibilidad magnética relativa y que es un número que no tiene dimensiones. Es 
constante si el medio es lineal. Si el medio es anisótropo pero lineal, M = Y H donde Y 
es ahora un tensor simétrico de orden dos. 
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También se observa este fenómeno en moléculas de sustancias diamagnéticas, como 
el benceno. 

Cuando en un material isótropo la susceptibilidad magnética varía de un punto a otro 
del medio, un campo magnético homogéneo produce una fuerza de traslación 


proporcional y paralela al gradiente espacial de esa susceptibilidad. 


Consideremos un volumen mesoscópico V, de un material biológico diamagnético. 
Como el tensor de susceptibilidad es simétrico, podemos orientar el sistema de ejes de 
coordenadas local según las direcciones principales de ese tensor. 

La matriz que lo representa resulta diagonal. 

Los elementos sobre la diagonal son las susceptibilidades principales Y, , X2 Y X3 , 


todas ellas números reales negativos. 


Xx 0 0 
En ese caso la matriz mencionada se puede escribir así: 0x2 0 (8.1) 
0. 0 2% 


El aporte a la densidad de energía libre de Gibbs local, debida a procesos magnéticos, 
cuando el campo crece a partir de cero a temperatura y presión constantes se puede 
expresar, en ausencia de fenómenos de histéresis, por 

— $, B-dH=-2u9 H? — po $, M-dH (8.2) 
El segundo término del miembro de la derecha de esta ecuación se puede interpretar 
como el aporte debido a la magnetización local del material. *' 


La magnetización se expresa como función lineal del campo H: M=XY:H 


A partir de esta aproximación se obtiene para la densidad de energía libre: 


> 
H> Es 


> 1 > a 
Im = Hof, M:dH=-5MH XX + H (8.3) 
Si las componentes del campo de inducción Ben dirección de los ejes principales del 


tensor son H, , H, y Hz resulta: 


1 
Ima ==> Moa Hf + X2 H5 Es de H3) (8.4) 


11 La energía libre de Gibbs para un sistema termodinámico magnetizado (en ausencia de 
fenómenos de histéresis) con energía interna U , entropía S, temperatura absoluta T, 


volumen V y presión p se define por G = U + pV-—TS-— f/f (18 : di) dV donde la 
integral triple se puede expresar como suma de dos términos: f[f 6 Ho H?) dV + 


Si (uo Le M- dH) dV La segunda integral se extiende sobre el sistema considerado, 
mientras que la primera se extiende al sistema y a su ambiente. La variación en la energía 
interna viene dada por: dU =T dS—p dV + fff (H : dB) dV (Guggenheim, 1967). 
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Supongamos, además, que las susceptibilidades principales toman dos valores: 


M=XM Y x=X=Xx, - Ambas susceptibilidades xy y x, son negativas, pero 
asumiremos que |x| > |x,|. 
Si H es la magnitud del campo externo, 0 es el ángulo que forma É con la primera 


dirección principal y y es el ángulo que forma la proyección H, de H (sobre el plano 

determinado por las otras dos) con la segunda dirección principal: H, =HcosO0 = Hy 
Hz = H send cosqy = H, cosqp Hz =H sen0 seng = H, senp 

La parte dependiente de la magnetización de la densidad de energía libre queda en este 


caso: 


Gui == Ho cos*0 + y, sen?0) B? ==> po H? Qu + Ga 24) cos?) 
El torque por unidad de volumen que actúa sobre el elemento de material biológico 
puede estimarse por la derivada parcial de gy y respecto del ángulo, a igualdad de las 
demás variables (temperatura, presión). 


Entonces el torque sobre la masa de estructura mesoscópica de volumen V, se puede 
calcular mediante la fórmula: Tm = => Hp H? V, Ax sen 20 (8.5) 


En esta última fórmula Ax =x.— 41 =1| Xx) —lx,1 > 0 

Ese torque por unidad de volumen tiende a rotar la estructura de modo que la dirección 
principal correspondiente a susceptibilidad de mayor valor absoluto se aproxime a la 
dirección del campo externo. 

Lo que acontezca en definitiva y por tanto los posibles efectos biológicos del campo 
homogéneo, depende de la relación de fuerzas entre las que producen el torque debido 
al diamagnetismo anisótropo y las demás fuerzas actuantes en la superficie y el volumen 
del elemento mesoscópico de material considerado. 

En estructuras con un grado elevado de anisotropía diamagnética y valores absolutos 
de los coeficientes de susceptibilidad también elevados, con un anclaje no muy fuerte a 
las estructuras vecinas, como por lo general es el caso de los cromosomas en el núcleo 
celular, se puede observar orientaciones inducidas por campos magnéticos 
homogéneos lo bastante intensos. Lo mismo se puede observar en membranas 
débilmente conectadas con las estructuras adyacentes. 

En general, cualquier molécula o estructura supramolecular con diamagnetismo 
anisótropo tenderá a orientarse de tal forma que la dirección de mínima susceptibilidad 
sea paralela al campo aplicado. El grado en que se pueda orientar un sistema 
diamagnético en un campo magnético externo a una cierta temperatura, en equilibrio 


termodinámico local, depende de una función de distribución angular que permite 
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estimar la probabilidad de que la dirección de y en la estructura forme un ángulo dado 
con el campo externo. 
La molécula de benceno presenta xy =—5.7 x 107? y x, =—0.15 x 107? siendo la 


dirección correspondiente a y, ortogonal al plano que se puede asignar a la molécula. 
En un campo externo la molécula experimenta un torque que tiende a alinear su plano 
con la dirección del campo. Para una molécula aislada que pueda rotar libremente en 
un campo externo, a temperaturas como las que se suelen encontrar en los sistemas 
biológicos, la agitación térmica se opone eficazmente al efecto debido al campo externo. 
Para un sistema compuesto por un gran número de moléculas, como las estructuras 
supramoleculares de las células o de los espacios extracelulares en los tejidos, el efecto 


de la agitación térmica disminuye significativamente. 


Pero aún en un medio que desde el punto de vista del magnetismo es isótropo, pueden 
aparecer fuerzas, en este caso de traslación, cuando la susceptibilidad magnética varía 
en el espacio. 

La fuerza magnética por unidad de volumen es proporcional al cuadrado del campo y al 


gradiente espacial de la susceptibilidad: f =- > HA?V xy (8.6) 


(9) Magneto-hemodinámica y resonancia iónica. 

La sangre es un medio conductor iónico que se mueve a través del sistema 
cardiovascular impulsada por la contracción del músculo cardíaco. En cada elemento 
de volumen mesoscópico se cumple el principio de electroneutralidad, de modo que hay 
un equilibrio local de carga entre los iones positivos y los ¡ones negativos, que se 
mueven formando parte de una solución acuosa que presenta componentes 
macromoleculares y celulares. 

En presencia de un campo magnético externo, una fuerza de Lorentz tiende a desviar a 
los iones en dirección perpendicular a su velocidad de movimiento conjunto con el fluido 


y perpendicular al campo. Si el ion posee una carga q y se mueve con velocidad % en 


un campo magnético de inducción B la fuerza de Lorentz sobre el ion F = q (5 eS B) 
tiende a separar las cargas de distinto signo que se encuentran en un mismo elemento 
de volumen y se mueven con la misma velocidad. 

El proceso de separación de cargas produce un campo eléctrico que se le opone. 
Consideremos un modelo simplificado de una arteria, que la asimila a un tubo cilíndrico 
de sección circular dentro del cual fluye la sangre. La sangre es considerada, a estos 


efectos, como una solución electrolítica en agua. 
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En estado estacionario aparece una diferencia de potencial V¿ en los extremos de un 
diámetro de longitud d perpendicular al campo externo, en una sección transversal al 
flujo: *? V, = vB d senó (9.1) 

En esta ecuación 0 es el ángulo que forman la velocidad representativa del flujo v y el 
campo magnético externo. 

Para el caso de la aorta, tomando v = 0.6 "/¿ ,d=0.025m y B=1T el valor 
máximo de la diferencia de potencial que se puede obtener es V¿ = 0.015 V. 

En general el diámetro arterial disminuye en dirección del flujo promedio de sangre, las 
arterias se ramifican y se curvan. 

La aorta presenta un tramo ascendente, un tramo curvo (el cayado) que se continúa por 
un tramo descendente (toráxico y abdominal) hasta que este último se bifurca para 
irrigar los miembros inferiores. En todo este trayecto se va ramificando. El flujo de 
sangre es pulsátil: tiene una componente constante combinada con una componente 
oscilatoria, de modo tal que la velocidad varía en una misma sección transversal durante 
el ciclo cardíaco. Además, la geometría de la aorta presenta variaciones de una persona 
a otra. 

Pese a que el modelo de arteria como tubo cilíndrico constituye una simplificación muy 
significativa, permite obtener algunos resultados que parecen correctos en orden de 
magnitud. 

Se han medido aumentos en las amplitudes de las ondas T del electrocardiograma 
cuando el organismo se encuentra en presencia de campos estáticos lo bastante 
intensos (campos con intensidades superiores a 0.3 T ya permiten medir ese efecto) y 
adecuadamente orientados. Hay cierta evidencia acerca de que este refuerzo en la onda 
T (que corresponde la repolarización del musculo ventricular y se acompaña de un flujo 
significativo en la primera porción de la aorta) se debe a la suma del potencial de 
repolarización con el potencial inducido por el campo magnético externo. 

La electroneutralidad a nivel mesoscópico, en ausencia de campos eléctricos 
adecuadamente orientados, elimina la posibilidad de una corriente eléctrica no nula en 
dirección del flujo y debida solamente al arrastre de los iones por el flujo de sangre. 
Las corrientes eléctricas aparecen perpendiculares a la dirección del flujo y no 
desaparecen hasta que se instale un campo eléctrico que equilibre el efecto del campo 
magnético. Esto último en general no ocurre debido a las variaciones que acompañan 


al flujo de sangre en el espacio y en el tiempo. 


12 Relacionado con el efecto Hall. 
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El marco teórico de la magnetohidrodinámica que en principio puede aplicarse en 
hemodinámica plantea una densidad de corriente eléctrica Jj= 0] (E +0 x B) en un 
punto de un fluido de conductividad eléctrica d , junto con la ley de Ampere V x H= J 
(despreciando la corriente de desplazamiento), la ley de inducción de Faraday V x É= 
-B y la ecuación del movimiento del fluido conductor de densidad ps , viscosidad 
dinámica y y presión p *? : ora v+ (5 15) = Vp+J]xB+nA0 (9.2) 
Si despreciamos los esfuerzos viscosos para simplificar, a partir de estas cuatro 


ecuaciones se obtienen dos nuevas relaciones: 


-Una nueva versión de la ecuación del movimiento del fluido conductor que involucra 


Lo se 1 ee D 
una presión magnética local rod que se suma a la presión del fluido y esfuerzos 
PE 1-73 => ; A AA 
mecánicos locales D (B .V) B- asociados a gradientes del campo magnético: 


er(8+ (6 -W)D)= -—v (p +82) +2 (B.v)B (9.3) 


-Una ecuación para la evolución del campo magnético: 

o 7 1 — ER —> 

ES RO x H) (9.4) 
La relación B= uH (con y = pp) y la aparición de la velocidad % del fluido en ambas 
ecuaciones las acoplan entre sí. Se introduce un coeficiente de difusión magnética 


Din — . A partir de un valor representativo de la velocidad del fluido v y de una 


longitud £ característica del sistema (por ejemplo, el diámetro arterial) se define un 


número sin dimensiones que se suele denominar número de Reynolds magnético!* 
EA ñ » . ñ ñ 

Rem = = . Cuando el orden de magnitud de este número es superior o inferior a 1, se 
m 


pueden hacer simplificaciones en las ecuaciones del modelo. 

Tanto los resultados de corridas de simulación digital basadas en modelos matemáticos 
numéricos como los de experimentos con animales o voluntarios humanos sugieren que, 
excepto en presencia de campos lo bastante intensos, con velocidades de flujo muy 
rápidas y en grandes arterias, no cabe esperar observar efectos magneto- 
hemodinámicos (Kyriakou et al., 2012; Martin et al., 2012). 


13 En esta y otras ecuaciones que involucran campos vectoriales, como por ejemplo un 
campo A (7) función de la posición T” : AA= V (v . A) =Vx (v Se A) 


14 Puede interpretarse también como una especie de número de Péclet como los que 
aparecen en la teoría de los procesos de transporte de masa por advección y difusión. 
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En campos magnéticos muy intensos (10 T) las fuerzas de arrastre debidas a efectos 
magneto-hemodinámicos pueden afectar ligeramente a la presión arterial y a la 


velocidad axial del flujo en una aorta de un modelo experimental animal (Markov, 2015). 


Si la fuerza F =q (3 x Ba) debida a un campo estático actuara sobre un ¡on no 


h . 7 . A ' mov 
restringido, de masa m , este se movería en una trayectoria circular de radio R = a 


Si el campo estático By se combinara con un campo alterno se produciría un estado 
resonante para una cierta frecuencia angular «wp del campo alterno: wW = - Bo 


(resonancia iónica de ciclotrón). 

Como la resonancia iónica de ciclotrón requiere una trayectoria no restringida bastante 
extensa, generalmente no disponible para el ion en una célula, no parece que este 
mecanismo pueda explicar la sensibilidad a la frecuencia del campo alterno que se 
observa de algunas reacciones bioquímicas que involucran al ion Ca**. Eso, pese a que 
las frecuencias de resonancia de esas reacciones coincidan, dentro del error 
experimental, con la frecuencia angular de ciclotrón del ion (218 Hz para By = 504 T). 
Lo que se plantea es un mecanismo diferente, denominado resonancia paramétrica 
iónica, que puede influir sobre los procesos en el interior de la célula a frecuencias 
angulares de resonancia paramétrica iónica wp que resultan ser, por mera coincidencia, 
cercanos a múltiplos de la frecuencia angular de resonancia iónica de ciclotrón (Polk y 
Postow, 2006). 

Pero a diferencia de esta última, dependen del cociente entre la amplitud del campo 
estático y la amplitud de oscilación del campo alterno. 

Generalmente las frecuencias angulares de resonancia paramétrica iónica se 


encuentran comprendidas entre 62.8 y 628 Hz. 


(10) Electrodinámica a nivel de tejidos biológicos y órganos. 


Cuando las relaciones constitutivas entre el campo de inducción B y el campo 
magnético, entre el campo de desplazamiento D y el campo eléctrico E y entre el campo 


de densidad de corriente eléctrica de conducción E y el campo eléctrico E, son lineales, 
las ecuaciones de Maxwell, a nivel tisular **, pueden escribirse de la siguiente forma: 


OB 


(10.1) VxE=-== (10.2) 


0D 


OE >» 0D 
VxH=)J Le ae 


15 Recordar la definición de los cuatro niveles en los que pueden ser estudiados los 
sistemas biológicos desde un punto de vista del bio-electromagnetismo, que aparece en 
el comienzo de la introducción al presente trabajo. 
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V:B=0 (10.3) V:D=p (10.4) 

B=ñ-H (10.5) D=é-E (10.6) J.=5-E (10.7) 
Aquí Ji es el campo de densidad de corriente eléctrica impuesta en el conductor de 
volumen debida a la actividad eléctrica de los tejidos (en general asociada a la génesis 
y a la propagación de potenciales de acción en las membranas de las células 
excitables), mientras que J. es la densidad de corriente pasiva, que aparece en 


respuesta al campo eléctrico local. El campo impuesto Ji , como lo sugiere su nombre, 
se considerará como un dato en el marco de las ecuaciones de Maxwell. 

Además, p es la densidad de carga libre en el punto considerado del conductor de 
volumen de los tejidos. 


poe > >. .0 ; : 
En la ley volumétrica de conservación de la carga V :] + =7 = 0 la densidad de corriente 


J' es la suma de la corriente impuesta con la corriente de conducción. 
Los campos se expresan en función del instante de tiempo t y en función del vector de 
posición Y . Como antes en este trabajo, el operador Vx representa el rotor del campo 


al que se aplica, y el operador V- representa su divergencia. 


La posible anisotropía del conductor de volumen se tiene en cuenta representando la 
permitividad dieléctrica, la permeabilidad magnética y la conductividad eléctrica por 
tensores simétricos, £ , il y O” respectivamente. 

En estos parámetros fenomenológicos el aporte de las células individuales no aparece 
en forma explícita: por este motivo la descripción se efectúa al cuarto nivel, el tisular y 
no en el tercer nivel, el nivel celular. 

A este nivel se consideran las membranas excitables, los espacios intracelulares y los 
espacios extracelulares como si fueran medios continuos superpuestos. Este tipo de 
modelos se conoce, en electrofisiología y biofísica, como modelos de bidominio. 
Cuando se consideran los campos magnéticos producidos por la totalidad de un órgano, 
como es el caso de la magneto-cardiografía o la magneto-encefalografía, la distribución 
espacial y temporal de las corrientes impuestas J! (+ ,7 ) en el conductor de volumen 
debida a la actividad eléctrica de los tejidos también se describe en el cuarto nivel, el 
tisular. 

En el caso del músculo cardíaco las células se encuentran acopladas eléctricamente en 
forma directa a través de conexiones que dan un cierto grado de continuidad al espacio 
intracelular. Esto, que se denomina algunas veces “sincicio eléctrico”, permite para 
diferentes propósitos (pero no en todos los casos) construir modelos útiles 


permaneciendo a nivel cuatro. 
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En el caso del sistema nervioso esa continuidad eléctrica por lo general no existe. Si 
bien la propagación de un potencial de acción a través del axón de una neurona puede 
influir en la polarización de las membranas de las neuronas vecinas, esa influencia por 
lo general no es decisiva. La conexión que posee la mayor importancia funcional es la 
que se establece a través de las sinapsis entre una célula nerviosa y otra. Es una 
conexión de índole electroquímica, que involucra la liberación y difusión de mediadores 
moleculares para conectar las células (hay sinapsis eléctricas, pero no son frecuentes). 
Esa ruptura de la continuidad eléctrica se puede tener en cuenta integrando modelos de 


nivel tres con modelos de nivel cuatro. 


> 


Los posibles efectos de histéresis en las relaciones constitutivas B = 1: H, 


D=¿:E y Jo = %- E, cuando los campos involucrados son oscilantes, se manifiestan 
como una variación en los parámetros fenomenológicos del modelo (fundamentalmente 
la permitividad dieléctrica y la conductividad eléctrica), con la frecuencia de oscilación. 
Los campos se descomponen en armónicas y se emplean expresiones complejas del 
tipo A(t, 7) = expliwt] A(P) paraE(t,F) , J(t,P) y H(t,P). 

Entonces los coeficientes de los tensores de permitividad dieléctrica, de conductividad 
y de permeabilidad magnética aparecen como funciones complejas de la frecuencia 


angular, en las relaciones constitutivas que relacionan los campos entre sí. 

La Figura 5 muestra la variación medida en la parte real de la permitividad dieléctrica 

relativa e, = - (eje de ordenadas a la izquierda) y la parte real de la conductividad 
(0) 


eléctrica (eje de ordenadas a la derecha) de un medio isótropo, en función de la 


frecuencia f = a Como es usual, ey representa la permitividad eléctrica del vacío. 


108 100 
o” 
e 104 
10 
10 
£ 
y 4 (09) 
10 p E 
4 
10? 
Y 
100 10? 101 108 108 1010 
f[Hz] 
Figura 5 
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Los resultados experimentales sobre esa variación con la frecuencia se interpretan 
mediante estructuras y procesos pertenecientes a los tres niveles jerárquicos previos 
(nivel atómico-molecular, nivel de las estructuras supramoleculares y nivel celular). 

Por debajo de los 100kHz las propiedades de conductividad de los tejidos están 
dominadas por la conducción de los electrolitos a través del espacio extracelular. La 
conductividad promedio del tejido depende fundamentalmente de la conductividad del 
medio extracelular y de la fracción del volumen del medio extracelular. 

Por encima de 100kHz la conductividad promedio depende también de la conducción 
de electrólitos a través del interior de las células debido a que desaparece 
progresivamente la barrera eléctrica establecida por las membranas celulares. 

Desde el punto de vista eléctrico las membranas biológicas se comportan como una 
combinación en paralelo de un capacitor con un resistor. Como la reactancia del 
capacitor de membrana tiende a cero a medida que aumenta la frecuencia de la 
corriente alterna, y el capacitor se halla en paralelo con la resistencia asociada a los 
canales iónicos de la membrana, la corriente alterna de alta frecuencia pasa a través 
del capacitor de membrana casi sin oposición, en forma de corriente de desplazamiento, 
eludiendo pasar a través de los canales ¡iónicos en forma de una corriente de 
conducción. Como consecuencia a altas frecuencias se produce un cortocircuito en el 


capacitor de membrana y se produce un aumento en la conductividad promedio. 


A bajas frecuencias el tejido muestra una dispersión dieléctrica que se llama dispersión 
a. y que se produce en el orden de los kHz. La dispersión a se debe fundamentalmente 
a dos procesos: a la polarización de contraiones próximos a superficies cargadas 
eléctricamente presentes en los tejidos y a la polarización de grandes estructuras 
macromoleculares que están fijas a las membranas. 

A frecuencias por debajo de las que corresponden a la dispersión a, la permeabilidad 
dieléctrica relativa alcanza valores de decenas de millones. La conductividad por su 
parte varía poco en la zona de frecuencias en la que se manifiesta la dispersión a. 
Entre 0,1 y 10 MHz, el tejido exhibe la dispersión f que se debe a la variación de la 
carga eléctrica en las membranas celulares que se produce por un desbalance entre la 
corriente eléctrica en el espacio intracelular y la corriente eléctrica en el espacio 
extracelular. 

Por encima de los 10 MHz la impedancia de las membranas celulares se hace 
despreciable y las corrientes pasan tanto por el espacio intracelular como el extracelular. 


Al pasar por la región correspondiente a la dispersión fP la permeabilidad dieléctrica 
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relativa disminuye significativamente pero también es notorio el incremento de la 
conductividad del tejido. 

A frecuencia de microondas por encima de 1GHz, el tejido exhibe una dispersión 
llamada y (que se centra en los 20 GHz). Se debe a efectos de relajación de la rotación 
de las moléculas de agua presentes en el tejido y es el mismo fenómeno que se 
encuentra en el agua líquida pura. Junto con la dispersión y en la permeabilidad 
dieléctrica relativa se observa un gran incremento en la conductividad del tejido. 

Tanto la conductividad como las propiedades dieléctricas de los tejidos dependen de la 


temperatura. 


Esta descripción de las dispersiones y las variaciones en la conductividad y la 
permeabilidad dieléctrica en un tejido blando es una descripción basada en valores 
promedio, tomados sobre el volumen ocupado por la muestra de tejido estudiada. 

Pero en general los tejidos biológicos comprenden materiales muy heterogéneos y en 
un organismo inmerso en un campo externo numerosos procesos ocurren en interfaces. 
Se encuentran grandes diferencias de conductividad entre, por ejemplo, los tejidos y 
fluidos de los vasos sanguíneos y el tejido conectivo que forma parte de los huesos. 
Aunque en menor grado, pero con importantes consecuencias sobre la estimulación 
funcional de fibras nerviosas, hay una diferencia de conductividad entre la sustancia gris 
y la sustancia blanca en el cerebro. Desde el punto de vista eléctrico entonces, no se 
puede considerar a los tejidos como un material homogéneo (Suárez-Antola, 2010 (a)). 


La conductividad eléctrica y la permeabilidad dieléctrica pueden depender de la 
dirección en la que se las mida, debido a la existencia de patrones morfológicos que 
introducen anisotropía, tales como haces de fibras musculares, tendones, haces de 
fibras nerviosas y vasos sanguíneos. Entonces, en menor o mayor grado todos los 
tejidos pueden presentar un comportamiento anisótropo que se describe mediante los 
tensores e y O. 

Para ilustrar la anisotropía en la conductividad eléctrica, se puede utilizar el esquema 
que aparece en la Figura 6. Es un modelo eléctrico muy idealizado de un tejido, cuyas 
células aparecen como elipsoides alargados con sus ejes mayores paralelos entre sí. 
Esto define una dirección preferencial. 

Se pueden ver dos pares de electrodos el par 1-2 que genera una corriente transversal 
a la dirección preferencial y el par 3-4 que genera una corriente paralela a dicha 


dirección. 
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Figura 6: Líneas de corriente asociadas 
a dos pares de electrodos en un modelo 
de tejido anisótropo. 


A bajas frecuencias, las líneas que conectan los electrodos sugieren las trayectorias de 
la corriente en este caso. La resistencia medida entre 1-2 es mayor que la resistencia 
medida entre 3-4 porque el camino que sigue la corriente transversal es más tortuoso. 
Como sugiere la Figura 7, además de ilustrar la anisotropía, este tipo de modelo 
idealizado ilustra uno de los mecanismos por los cuales la conductividad eléctrica puede 


variar con la frecuencia de la corriente alterna. 


ze o Aé e e 


Figura 7. Caminos de la corriente en un 
tejido para baja (LF) y alta (HF) 
frecuencia. 


HE. connnneas LE — 


A frecuencias lo bastante elevadas la reactancia de las membranas celulares y de las 
membranas de los orgánulos intracelulares prácticamente desaparece debido a que, 


como dijimos previamente, se cortocircuita el capacitor de membrana. 


Desde el punto de vista magnético, el comportamiento de los tejidos biológicos, 
considerados globalmente, es por lo general mucho más simple que su comportamiento 
desde el punto de vista eléctrico. Si Í representa al tensor identidad y F representa al 
tensor de susceptibilidad magnética: Hi = o (+) (10.8) 

Como las componentes del tensor Y en los tejidos biológicos en masa, diamagnéticos, 
son, en valor absoluto, varios órdenes de magnitud numérica inferiores a 1 (del orden 


de 10* o menor) se puede utilizar la simplificación fi = uy Y que implica considerar los 
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tejidos como magnéticamente isótropos y desde el punto de vista magnético, 
equivalentes al medio el cual el organismo correspondiente se encuentra inmerso. 
Los magnetosomas, en caso de ser necesario tenerlos en cuenta, deben tratarse por 


separado como estructuras mesoscópicas, es decir, a nivel dos. 


En los modelos electrodinámicos a nivel de tejidos biológicos se emplean las conocidas 
condiciones en la frontera para los campos y las densidades superficiales de carga y de 
corriente de la teoría electromagnética. 
En una interfaz entre dos medios o dos fases de un tejido considerado como un medio 
polifásico, digamos la fase 1 y la fase 2, la componente normal de la inducción magnética 
y la componente tangencial del campo eléctrico son continuas: 

Bin = Ban (10.9) Et = Eze (10.10) 
La diferencia entre las componentes normales de la inducción eléctrica a uno y otro lado 
de la interfaz es igual a la densidad superficial de carga eléctrica w, que pudiera hallarse 
en la interfaz: Din — Dan = 0s (10.11) 
Si la densidad de corriente eléctrica en la interfaz es nula, la componente tangencial del 
campo magnético es continua en la interfaz: Hi = Has (10.12) 
La diferencia entre las componentes normales de la corriente eléctrica a través de la 
interfaz verifica una versión de superficie de la conservación de la carga: 

Jan Jan + 005 =0 (10.11) 

Las regiones abarcadas por las corrientes inducidas incluyen varios tipos de tejidos 
diferentes e interfaces entre unos y otros. Por lo general no resulta factible tener en 
cuenta los detalles de la distribución espacial de las propiedades eléctricas del tejido en 
masa y de sus interfaces. 
Aún con el apoyo de modelos numéricos, es necesario trabajar con valores 


promediados. 


(11) Inducción de corrientes eléctricas en el conductor de volumen de 
los tejidos: disipación de energía y activación de tejidos excitables. 


Todos los sistemas de cargas aceleradas radian energía electromagnética, es decir 
actúan como antenas emisoras, pero no todos esos sistemas lo hacen con la misma 
eficiencia. 


Si 1 es una longitud de onda característica emitida por un sistema de cargas aceleradas, 
se pueden definir dos regiones tomando como referencia la longitud alo 7, * Cuando la 


distancias entre un punto de observación y el sistema emisor, es sensiblemente inferior 
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a al 7 » domina el patrón de campo electromagnético denominado campo cercano, 


mientras que si esa distancia es sensiblemente superior a o domina el patrón 


denominado campo lejano o campo de radiación. 

El campo cercano intercambia energía electromagnética con la fuente del campo, 
mientras que el campo lejano propaga como una radiación la energía electromagnética 
a partir de la fuente del campo, sin retornarla, a menos que sea convenientemente 


reflejada por algún objeto lejano. 


A 50 Hz la longitud ada ,, en el aire es de 960 km. A 100 kHz, de es de 4800 m. A1 


MHz, de 7, toma el valor 48 m. A centenares de GHz, de 7 Vale fracciones de milímetro. 


Un organismo puede interactuar y extraer energía tanto si se encuentra en el campo 
cercano como en el campo lejano, pero los mecanismos de interacción son diferentes 
en uno y otro caso. 

Supongamos que un organismo se encuentra en el aire, en el campo lejano, a una 
distancia grande respecto de la frontera entre el campo lejano y el campo cercano de 
una fuente de ondas electromagnéticas. 

En ese caso, sobre el organismo incide una onda prácticamente plana donde el campo 
eléctrico y el campo magnético son ortogonales entre sí, oscilan al unísono, 


indisolublemente vinculados, uno con otro. Sus magnitudes se relacionan a través de la 
fórmula: ==3770 (11.1) 


La penetración del campo electromagnético desde la superficie hacia el interior del 


organismo se produce como un fenómeno de propagación de una onda que involucra a 


0D ] a : 
la corriente de desplazamiento se y Ya asociada con una disipación de energía en 


forma térmica que puede ser muy significativa. 

Como los efectos térmicos se deben a las corrientes inducidas, y no a los campos 
eléctricos o magnéticos interactuando directamente con las estructuras tisulares, en 
principio esos efectos pueden ser producidos por campos de frecuencias muy 
diferentes. 

La densidad de energía disipada en un punto de un conductor de volumen, que se 
convierte localmente en energía calorífica viene, dada en función de la densidad de 


corriente eléctrica inducida o del campo eléctrico local, por: 
=P =0E (11.2) 
La Figura 8 muestra los resultados de una simulación digital de la distribución de la tasa 


de absorción específica dentro de la cabeza de un joven irradiado por un teléfono celular 
a 902 MHz. 
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Figura 8. Resultados numéricos debidos a Myoung Soo Kwon (Lin, 2012) 


La tasa de absorción específica aumenta a medida que ascendemos en la barra vertical 
pasando del azul oscuro hasta el blanco en el extremo superior (allí la tasa de absorción 
específica es de 1.25 W/kg). A la izquierda de la cabeza que aparece en la figura, a la 
altura del oído izquierdo, se puede apreciar la ubicación del teléfono celular y su antena. 


Para un organismo, si la totalidad de la energía calorífica permanece en el punto en el 
que es entregada, la velocidad de variación de la densidad local de energía interna se 
relaciona con la velocidad de variación de la temperatura local T , con la densidad local 


de corriente eléctrica J y con la conductividad local 4 mediante la expresión: 


0 A 
Pm Cr = 734 (Kp VT) +F, ps *C5 (Ts —T) + Gmet (11.3) 


Suponiendo que se conocen los valores locales de la densidad másica p,,, del calor 
específico c , de la conductividad eléctrica y de la conductividad térmica K+y local del 
tejido, el flujo de perfusión local de sangre F; , la densidad p, , el calor específico cs 
de la sangre, su temperatura T, y la densidad local de generación de calor por el 
metabolismo qme+ , esta ecuación junto con las condiciones necesarias iniciales y en la 


frontera, posibilita estimar el incremento de temperatura que se produce debido a la 


conversión en calor de la energía asociada a la corriente eléctrica. 
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Si los parámetros se suponen constantes y se desprecia el transporte de calor desde el 
lugar en el que aparece hacia las regiones vecinas y hacia la sangre, así como el aporte 
debido al metabolismo, se tiene para el incremento AT de temperatura luego de un 


P+dmee)to (114) 


Esta expresión puede utilizarse para estimar una cota superior al incremento de la 


intervalo de tiempo de breve duración t: AT = E a 


temperatura en un punto del conductor de volumen del organismo, puesto que el flujo 
de calor hacia otros lugares, por conducción o por convección a través de la sangre, 
tiene como consecuencia disminuir AT en igualdad de las demás condiciones. 

El efecto de oír microondas, descubierto durante la segunda guerra mundial, resultó ser 
en última instancia de origen térmico. Consiste en sensaciones auditivas que una 
persona experimenta cuando su cabeza está expuesta a microondas pulsadas, como 
las generadas por un radar. Un muy pequeño incremento de temperatura que se 
produce muy rápido en el tejido cerebral genera por expansión termo-elástica una onda 
acústica que viaja hacia la envoltura craneana, se refleja y transmite, y finalmente llega 
al oído interno siguiendo un proceso auditivo normal. 

La intensidad percibida del efecto acústico, para una potencia incidente constante, 
depende de la duración del pulso de radar, observándose incrementos y disminuciones 
en la intensidad percibida asociadas a resonancias que se producen en la cavidad 
acústica equivalente a la cabeza humana (Reilly, 1998). 

Actualmente se considera que el principal mecanismo conocido por el cual las emisiones 
de radio y microondas producen efectos biológicos es el calentamiento de los tejidos. 
En experimentos con animales se observan daños en los tejidos y otros efectos 
inducidos térmicamente cuando la cantidad de energía absorbida por el animal excede 
significativamente la cantidad de calor generada por los procesos corporales normales. 
Los efectos térmicos de las radiofrecuencias o de las microondas, al mismo nivel de 
potencia que la generación de calor corporal en condiciones basales, parecen no 
representan riesgos biológicos. 

Esto no descarta la posibilidad de que existan efectos no térmicos, que podría aparecer 
como resonancias para ciertas ventanas de frecuencia. 

Para determinar en qué medida en un organismo biológico se podría producir una 
resonancia, el concepto de dimensión eléctrica de un sistema material resulta mucho 
más útil que las dimensiones geométricas tomadas por separado. 

La dimensión eléctrica es un número sin dimensiones que se define en relación con una 
oscilación de una frecuencia dada f en el material del cuerpo en el cual la velocidad de 


propagación de las ondas a esa frecuencia esc. 
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Se debe tener en cuenta que la frecuencia es la de la oscilación del campo externo, pero 
la velocidad de propagación es la del material del cuerpo. 
Si f es una dimensión espacial característica del cuerpo, y si 4 es la longitud de onda 


de la oscilación provocada por el campo externo en el material, la dimensión eléctrica 


se define así: d(f) = - = - (11.5) 


Para un cuerpo con una dimensión eléctrica adecuada en presencia de un campo 
electromagnético externo, se puede producir una resonancia electromagnética: en un 
cierto lugar del conductor de volumen formado por los tejidos se puede instalar un 
campo eléctrico local significativamente mayor que el que aparece en ese mismo sitio 
para frecuencias de oscilación ubicadas fuera de la resonancia. 

Como en condiciones de resonancia las estructuras biológicas que se encuentran allí 
localizadas interactúan con un campo oscilante de mayor amplitud, parece más 
verosímil esperar una transferencia de energía con consecuencias sobre los procesos 
metabólicos y la transferencia de señales tisulares locales y eventualmente, globales, 


aunque los efectos térmicos no sean significativos. 


Si un organismo se encuentra dentro del campo cercano de un sistema emisor, lejos de 
su frontera con el campo lejano, los efectos debidos a la propagación del campo 
electromagnético en general pueden ser despreciados: la corriente de desplazamiento 


0D ST 
2 eE puede despreciar. Entonces, en el interior del organismo la ley de Ampere- 


Maxwell se reduce a: VxH=]Ji+]. (11.6) 
Las demás ecuaciones de la electrodinámica se mantienen incambiadas, excepto la ley 
volumétrica de conservación de la carga que ahora se reduce a V -] =0. 

En este caso la disipación de energía en forma térmica no suele ser significativa y los 
efectos de los campos externos que se pueden observar se relacionan por lo general 


con la excitación de nervios, músculo esquelético y miocardio. 


Teniendo en cuenta los dispositivos utilizados en el tratamiento y en el diagnóstico 
médico, y las características de los campos producidos por la actividad fisiológica de los 
organismos, consideraremos de ahora en adelante frecuencias comprendidas entre O y 
100 kHz. 

Consideremos un modelo muy simple de dos medios, aire (medio 1 de conductividad o, 
y permitividad dieléctrica e, ) y un conductor de volumen isótropo (medio 2 de 


conductividad o, y permitividad dieléctrica e, ). 
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En la interfaz E¡¿ =E>¿ , mientras que Din — Dan = €1 Ein — €2Ezn =Ws Y 
0 ; 
Jin Tan = 01 E1n — 0 En == 7 %s Suponiendo que se ha alcanzado un estado 


estacionario de oscilación a una frecuencia w, un factor expli wt] aparece en los 


campos en ambos medios y en la densidad superficial de carga. 
0 : 8% 
Entonces os = 1005. Introduciendo esta relación en las dos ecuaciones para la 


componente normal del campo a uno y otro lado de la interfaz: 


¿wE1E1n — 10 €2E7n =100s 01E¡n — 07E23n =-100s 
De estas ecuaciones resulta: (0, + ¡we Ejn = (0, + ¡wez) E>n 
. ; E 
Si y es el ángulo que el campo forma con la normal a la interfaz, tenemos tg y = = 
n 


Entonces, teniendo en cuenta la relación entre E; , y E» n obtenemos para los ángulos 


an) tg pa (11.7) 


a uno y otro lado de la interfaz: tg Pp = (E na 
Para frecuencias inferiores a los MHz, esto implica que el campo E, en el aire tiene casi 
la dirección de la normal a la interfaz: el organismo distorsiona el campo eléctrico en el 
aire de modo que incida casi perpendicularmente a la interfaz. ** 

La Figura 9 esboza esta situación para un campo eléctrico oscilante que es vertical y 


homogéneo lejos de la interfaz con el cuerpo humano. 


16 Para el aire o, = 1078 $/,, e, = 1071 F/, mientras que para el conductor de volumen 
de los tejidos vamos a tomar od, = 1 DS EA y e, =10* ce para una frecuencia angular 


w = 10% Hz. Resulta arras = 107” de modo que tg p, = 107 tg p, Entonces, si tg p, = 
2 2 


100 (un campo E, en el conductor de volumen, que forma un ángulo apenas menor que 
T/, con la normal a la interfaz) todavía se tiene tg p, = 10”. 
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Si 


7. 


Figura 9. Distorsión, debida al cuerpo humano, en las líneas de un campo eléctrico 
oscilante en el aire. El sentido de las líneas se invierte cada medio ciclo de oscilación. 
Las líneas de campo interno y densidad de corriente eléctrica en el conductor de 
volumen no están a escala. (Modificado de Reilly, 1998). 


Como la longitud de onda a frecuencias inferiores al MHz es muy grande respecto de la 
máxima distancia entre puntos del cuerpo humano, se puede suponer un campo externo 
variable en el tiempo, pero constante en dirección y amplitud en todas partes del espacio 
próximo al organismo *”, siendo w = 2rr. f la frecuencia angular: 

Eaire(t) = Eycos[(w.t) (11.8) 
Considerando para simplificar una losa transversal al campo, de conductividad d(w) y 
constante dieléctrica relativa e, (w) , que se comportan en función de la frecuencia 
angular de la misma manera que se comportan esos mismos parámetros en un tejido 
biológico tipo. Supondremos que la losa se puede considerar infinita en sentido 
perpendicular al campo oscilante. 
Para los fines del presente modelo, se puede suponer que el aire que rodea la losa 
posee conductividad y = 0 y constante dieléctrica relativa e, = 1 
La simetría de la situación física permite suponer que el campo en la losa es paralelo al 
campo externo, o sea es perpendicular a las caras planas de la losa, como sugiere la 
figura 9. 


17 Los resultados de aplicar modelos más realistas (Lin, 2012; Polk y Postow, 2006, 
volumen 1; Greenebaum y Barnes, 2019, volumen 1) no cambian las conclusiones más 
importantes que se obtienen del modelo simplificado. 
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campo en el aire 


A o Í 
campo en los tejidos 
HB h 


Figura 9 


Si Etejidos(t) representa el campo en el interior de la losa que representa en este modelo 
alos tejidos y si ws(t) es la carga eléctrica inducida en una de las dos interfaces losa- 
aire (la carga inducida en la otra interfaz es —ws(t) , de la condición de borde para la 
componente del campo eléctrico normal a la interface se obtiene: 

tE aire (1) + Eo€rErejidos(t) = ws(t) (11.9) 
La ley de conservación de la carga eléctrica suministra la relación siguiente (teniendo 


en cuenta que solo hay corriente de conducción en el interior de la losa, pero no en el 


. dws 
aire): OErejidos = ==e (11.10) 


De estas dos ecuaciones se desprende: 


dEtejidos 


E + Etejiaos = — E Eo Sen(ot) (11.11) 


Er€o 


Luego de un transitorio cuya duración es del orden del tiempo de relajación 7, = 
se alcanza un estado estacionario de pulsación del campo eléctrico en el interior de la 
losa que representa a los tejidos biológicos. 

Estimando el tiempo de relajación para y = 0.5 Sm y €, = 50 resulta 1, =107?s. 


El campo oscilante estacionario  Etojigos(t) = Acos(wt) + Bsen(wt) se instala 


WT 


entonces muy rápido. Los coeficientes son: A = —(wT¿)B B= oda Eo 


Como para frecuencias hasta del orden de los MHz, wr, es mucho menor que uno, se 
puede despreciar A frente a B y se puede aproximar B por — Es a 


A partir de esta aproximación se obtiene que el campo oscilante en el interior de la losa 


viene dado por  Etejidos(t) Y — Es sen(wt) 
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De este último resultado se desprende que la densidad de corriente de conducción 


Jejidos = 9 * Etejiaos Pulsa con frecuencia «wW y posee una amplitud de oscilación 


we7yE,y independiente de la conductividad de la losa. ** 

Para una frecuencia de 50 Hz el campo E, es de un orden 10% veces menor que el 
campo E, en el aire. Para una frecuencia de 1 MHz la amplitud del campo disminuye a 
menos de la centésima parte. 

Si el campo externo es estático, es decir si w=0, el campo eléctrico en el interior de la 


losa es nulo, como podía esperarse en el caso de un cuerpo conductor en equilibrio. 


A frecuencias menores a 100 kHz la membrana celular apantalla al citoplasma respecto 
de los campos eléctricos que puedan presentarse en el espacio extracelular. Para una 
frecuencia de 50 Hz el campo justo afuera de la célula es 10% veces menor que el campo 
en el aire, en la membrana celular es 10? veces menor y en el interior celular es 10* 
veces menor (Lin, 2012; Polk y Postow, 2006; Greenebaum y Barnes, 2019). 

En principio, no cabe esperar que a este nivel el campo eléctrico en el citoplasma pueda 
afectar los orgánulos funcionales de la célula o al ADN. Por otra parte, se sabe ahora 
que tanto corrientes eléctricas débiles inducidas por un campo magnético generado por 
una bobina como por un campo eléctrico generado por un sistema de electrodos pueden 
contribuir a soldar una fractura de un hueso (Polk y Postow, 2006; Markov, 2015). 


No obstante, este efecto de pantalla se debilita a medida que crece la frecuencia. A 
frecuencias mayores a 1 MHz desaparece, porque la membrana es cortocircuitada a 
nivel de su capacitor. El único efecto de pantalla remanente es el que pueda producir el 
conductor de volumen de los tejidos considerado en bloque. 

A frecuencias inferiores a los 100 kHz, el campo eléctrico externo no penetra casi en los 
tejidos biológicos, por lo cual los efectos observados se deben fundamentalmente al 
campo magnético externo, que penetra sin dificultad debido a que, como vimos más 
arriba, las permeabilidades magnéticas del aire y de los tejidos toman valores muy 
próximos. 

La Figura 10 sugiere lo que ocurre cuando un organismo se encuentra inmerso en un 


campo magnético externo oscilante. 


18 Si este resultado se intentara describir desde el punto de vista de la teoría de circuitos 
eléctricos, habría que asumir que la fuente de fuerza electromotriz (debida al campo 
eléctrico externo) posee una impedancia interna tan elevada que impone una corriente 
independiente de la resistencia eléctrica (el conductor de volumen de los tejidos, 
representado por la losa) a la que alimenta. 
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TZ ZA 


Figura 10. En ambas figuras se esbozan las líneas de campo eléctrico en el interior del 
cuerpo humano inducidas por un campo magnético externo oscilante y homogéneo. A la 
izquierda el campo magnético es vertical y a la derecha es horizontal, como viene 


indicado por las flechas gruesas. (Modificado de Reilly, 1998). 


El campo eléctrico inducido se conecta con el campo de inducción magnética a través 


de la versión integral de ley de Faraday: 

SE(,7) 55) ds =- JH B(,F) -5(P) as (11.12) 
En la parte izquierda de la Figura 10 el campo B se puede considerar como un campo 
variable con el tiempo, vertical y homogéneo: B(t,F)=-—Bo(0) é, (11.13) 
Suponiendo que $E (t, 7) .£(P) ds se toma sobre una circunferencia horizontal de 
radio r incluida en el conductor de volumen y que el campo eléctrico es tangencial y de 
igual magnitu E(t, 1) = E(t, 7) .2(F) en cada punto de la circunferencia en un mismo 
instante de tiempo: $E (t, 7) .t(7) ds = 27 r E(t,r) (11.14) 
Por otra parte 7 (F) = é, y - B(t, 7) = -EBo (t) é, con lo cual la integral de 
superficie tomada sobre el círculo de radio r delimitado por la circunferencia 


mencionada es: z E B(t,F) ¡(PH dS =nr? LB, (0) (11.15) 
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Igualando ambas integrales se obtiene: E(t,r) = > r —B, (t) (11.16) 


Si se emplea un modelo de conductor lineal, homogéneo e isótropo de 


conductividad gd, la densidad de corriente inducida también es tangencial a la 


circunferencia y su magnitud viene dada por: — J(t,r) = > gr —B, (E) (11.17) 


Las fórmulas obtenidas para las magnitudes del campo eléctrico y de la densidad de 


/ l / a ' 
corriente sugieren que para una misma ley de variación Bo (t) esas magnitudes 


aumentan como el radio de la circunferencia. Por ende, este modelo sencillo sugiere 
que, durante una exposición de cuerpo entero a un campo externo homogéneo 
vertical, los máximos de campo eléctrico y de densidad de corriente se encontrarán 
sobre el círculo de mayor radio que se pueda dibujar, perpendicular al vector de 
inducción magnética, completamente incluido en el conductor de volumen formado por 
el cuerpo. 

Por el contrario, cuando el campo es horizontal y se ubica como se sugiere en la parte 
derecha de la Figura 10, el resultado de una exposición de cuerpo entero no se puede 
describir mediante un modelo sencillo. 

En general los modelos puramente analíticos más o menos simplificados (como el que 
representa el cuerpo como si fuera un elipsoide prolato con eje mayor vertical) no son 
aplicables para estudiar los campos inducidos por exposiciones de cuerpo entero a 
campos magnéticos homogéneos, como las que se producen durante los estudios por 
resonancia magnética. Es necesario recurrir a modelos numéricos lo suficientemente 
detallados (Lin,2012). 


La situación es muy diferente durante la estimulación funcional destinada a tratamiento 
médico o a diagnóstico. La Figura 11 ilustra este aspecto de la estimulación funcional. 


Magnetic 
field B 


———  Stimulating 
Coil 


Electric a 


volume 
field E conductor 


Figura 11. Modificada de la tesis de 
doctorado de Jarmo Ruohonen (1998) 
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Ahora se crean campos magnéticos localizados, mediante bobinas que se sitúan 
próximas a la zona que se desea estimular. Ya no son campos que se puedan considerar 


homogéneos respecto de las dimensiones del cuerpo, pero sí son casi-estacionarios. 
Se pueden expresar como el producto de un campo vectorial función de la posición F(F) 
por una función escalar f (t) del tiempo: 

BF) =f(0 E(7) (11.18) 
La función f(t) se puede identificar con la corriente i, (t) que se hace circular a través 


de la bobina estimuladora. 
El campo eléctrico local inducido viene dado por: 

É(t,P) = E le (9) 46) (11.19) 
Las corrientes inducidas generan a su vez un campo magnético propio, pero como se 
verá en la última parte de este trabajo, a frecuencias inferiores a 100 kHz esos campos 
secundarios pueden ser despreciados. 
Ahora es necesario considerar como se relaciona este campo eléctrico con la 
polarización de las membranas excitables de células nerviosas y musculares hasta 
llegar a generar la activación de la célula, manifestada por la emergencia de un potencial 
de acción capaz de propagarse. 
Para hacer esto, partiremos de un modelo matemático simple de una fibra excitable, 
asimilada a un cable delgado, posiblemente curvado, cuyo interior corresponde al 
espacio intracelular. La membrana celular lo separa de su ambiente, formado por el 
resto de los tejidos. 
A partir de un punto de la fibra, tomado como origen, introducimos una coordenada de 
posición s . Esta coordenada se aplica a una fibra curvada y corresponde a la longitud 
del arco de curva (orientado) que conecta el origen con el punto considerado. 
La conexión entre el campo eléctrico inducido por el campo magnético variable y la 
polarización de la membrana excitable se establece a través de la función activante. 
Esta es una función, definida para cada instante de tiempo y para cada valor de la 
coordenada de posición, como la derivada respecto de esa coordenada de la 


componente del campo eléctrico tangente a la fibra E(t,F (s)) - £(s) : 

O (2. > > d. dl(37r> > 

2 (E(6,7 (5) (59) = (E ip(0)) L (E (* (5) (5) (11.20) 
El vector unitario £(s) es tangente en cada punto s de la fibra. Si esta última cambia de 


dirección en el espacio, Es) +0 y aunque F (5 (s)) sea constante a lo largo de la 


fibra, como podría ocurrir en ciertos casos durante una exposición de cuerpo entero, la 
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función activante no va a ser nula y va a producir una polarización en la membrana de 
la fibra (Suárez-Antola, 2010 (b); 2012; 2013). 

La generación de un potencial de acción en la membrana de la fibra se produce cuando 
en esa membrana se alcanza un patrón de polarización umbral. 

Para ello una longitud mínima de membrana de una fibra, que podemos suponer se 
encuentra inicialmente en reposo, debe ser despolarizada por encima del denominado 
umbral de membrana. Esa despolarización resulta de corrientes iónicas que cruzan la 
membrana desde el interior de la fibra. 

Pero como no hay fuentes de corriente en el citoplasma que puedan originar este tipo 
de corrientes despolarizantes, la conservación de la carga eléctrica exige que las 
corrientes despolarizantes sean compensadas por corrientes iónicas que atraviesan la 
membrana, en otras partes, hacia el interior de la fibra. 

El estado umbral se alcanza o no dependiendo de cómo se genere ese patrón de 
corrientes (Suárez-Antola, 2010(b); 2012; 2013). 

Las características del patrón generado dependen de cuánto dure y de la amplitud del 
estímulo, así como del estado inicial de la membrana excitable. La eficacia del estímulo 


depende de las características del campo inducido (E lp (0) E (7) manifestadas a 


z .z . .z 0 E > 2 
través de la función de activación [E(8,7 (5) -t(9). 


La Figura 12 muestra un pulso de corriente generado en una bobina estimuladora de un 
equipo de estimulación magnética funcional de tercera generación junto con un patrón 
de magnitud de campo eléctrico inducido correspondiente a ese pulso (Peterchev et al., 
2008). 
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Figura 12. (Tomado de Peterchev et alt., 2008). 
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Un pulso como el que muestra la Figura 12 es un pulso umbral cuando el campo que 
aparece abajo tiene una amplitud y una duración suficientes. A menor duración, mayor 
tiene que ser su amplitud para alcanzar o superar el umbral. 

Los resultados obtenidos se pueden graficar en un plano, ubicando la duración en las 
abscisas y la intensidad en las ordenadas. Los pares duración-amplitud umbral forman 
una curva que se asemeja a una hipérbola con una asíntota vertical en el eje de 
ordenadas y una asíntota horizontal, paralela al eje de abscisas, que corresponde a una 
amplitud mínima, como sugiere la Figura 13. Si esa amplitud mínima, conocida como 
reobase, no es superada, la excitación fracasa por larga que sea la duración del pulso 
estimulador. La duración de un pulso cuya amplitud umbral es el doble de la reobase se 


denomina cronaxia. 
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Figura 13. 


La curva umbral separa en primer cuadrante en dos regiones: una región de estímulos 
que no alcanzan el umbral, debajo de la curva, y una región de estímulos que superan 
el umbral, encima de la curva. 

La determinación del umbral de estimulación debe hacerse en cada caso, porque la 
curva intensidad-duración depende de la posición y tamaño de las bobinas, de las 
propiedades del conductor de volumen formado por los tejidos biológicos y de la 
naturaleza y posición de la o las células blanco de la estimulación. 

Se tienen en cuenta pequeñas variaciones aleatorias introduciendo un criterio 
estadístico para determinar el umbral. Por ejemplo, se considera pulso umbral cuando 
la excitación resulta exitosa en el 50 % de los casos para una misma amplitud y duración 


prefijadas en el equipo estimulador. 


(12) Medición de los campos magnéticos generados por los procesos 
fisiológicos. Sus fundamentos en la teoría electromagnética de los medios 
continuos: aplicación del teorema de reciprocidad. 
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La actividad fisiológica del sistema nervioso (central y periférico), de los músculos esqueléticos 
y del corazón se acompañan de corrientes eléctricas impuestas JE (t,F) al conductor de 
volumen formado por los tejidos. Estas corrientes generan campos magnéticos tanto en ese 
conductor de volumen como en la región que rodea al cuerpo humano. 

Esas corrientes y campos varían con el paso del tiempo. 


Ubicando en un sitio apropiado las bobinas de un equipo receptor de señales magnéticas, las 
variaciones en los flujos de inducción que atraviesan las bobinas inducen en ellas fuerzas 
electromotrices que se pueden medir. 


Como las relaciones señal-ruido son muy desfavorables, en el caso de los fenómenos magnéticos 
asociados a los procesos fisiológicos se emplean equipos superconductores (SQUID) y 
mecanismos para eliminar ruidos*? con el fin de poder registrar magneto-cardiogramas, 
magneto-encefalogramas y la magneto-miografías. 


Veamos ahora un teorema de reciprocidad demostrado en 1895 por Hendrik Lorentz (Papas, 
1988), a partir del cual Robert Plonsey (Plonsey,1972) pudo relacionar la fuerza electromotriz 
inducida en la bobina receptora con el campo de las densidades de corrientes generadas en los 
tejidos por la actividad fisiológica y con el denominado campo de derivación. 

Este último campo se obtiene imponiendo una corriente eléctrica variable a los devanados de la 
bobina receptora, de forma tal que la variación temporal del campo magnético generado por la 
bobina induzca un campo eléctrico en los tejidos. 

El campo de derivación se define a partir de este campo eléctrico inducido en los tejidos, como 
veremos a continuación. 


Comencemos replanteando las ecuaciones de Maxwell, a nivel tisular. Como vimos en la sección 
sobre electrodinámica a nivel de tejidos biológicos y órganos, esas ecuaciones pueden escribirse 
de la siguiente forma: 


Vx H=P+j + vxE=-2% v-B=0  V-D=p 


Las relaciones constitutivas se pueden asumir lineales, aunque en general hay que tener en 

cuenta cierto grado de anisotropía y heterogeneidad, expresado por tensores simétricos de 

permeabilidad magnética fi , permitividad dieléctrica £ y conductividad eléctrica % : 
B=ñi:H D=%-E J,=3:-E 

La heterogeneidad puede describirse mediante campos tensoriales dependientes del punto del 

espacio Y considerado o introduciendo un número apropiado de fases homogéneas. 


$ 
Supongamos que en la región R, un campo de densidades de corriente impuestas / | genera en 
>» > 
el espacio los campos E, y H; . 
ez Ñ . A Pi á 
En otra región R¿ un campo de densidades de corriente impuestas J ; genera en el espacio los 


campos E, y H, . 


19 Blindajes, gradiómetros o conjuntos distribuidos de bobinas conectadas a sistemas de detección multicanales. 
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Supongamos, además, que ambos campos de densidades de corriente oscilan en forma 
monocromática a la misma frecuencia angular w : 


J ¡(t,F) = expliwt]Ji(w;7) J5(t,F) = expliwt]J (w;r) 


En estado de oscilación estacionaria, los campos eléctrico y magnético oscilan a la misma 


frecuencia: 
E, (t,F) = expli w t] E,(w; 7) H,(t,7) = expli w t] H,(w; F) 
E,(t,F) = expli w t] E,(w;7) Ha(t,F) = expli «w t] H¿(w; 7) 


Ahora todos estos campos vectoriales poseen componentes complejas: son denominados 
fasores vectoriales. 


Para un medio anisótropo (con tensores de anisotropía simétricos) y heterogéneo, el teorema 
de reciprocidad de Lorentz permite concluir la siguiente igualdad entre integrales de volumen: 


Sy Elo; 7) + Ji(0;7) aV = f, Er(o;F) - J¿(o;7) av 


El producto que aparece en ambos integrandos es el producto escalar de vectores y ambas 
integrales se extienden a todo el espacio físico, tomando como sistema de referencia el 
consultorio o el laboratorio y empleando como modelo el espacio euclídeo tridimensional R . 

Como Filo; 7) es nulo excepto en R; yJi(w; Y) es nulo excepto en R, , se obtiene la versión 
del teorema de reciprocidad que utilizaremos en lo que sigue: 


pad > Fi A? v= E? > E Fi > 
ds E2(w; Y): Ji(w; r) dV a Sa, E, (w; f) J2(w;r) dv 
Esta ecuación se verifica en una situación muy general, incluso para los casos en los que se radian 
ondas electromagnéticas desde las regiones R; y R,. 
Se puede construir un sistema bipolar de registro y estimulación magnéticos análogo a un 
sistema de registro y estimulación eléctrico basado en dos electrodos. 


Siguiendo una sugerencia de Baule y Mc Fee (1965) en el caso magnético las superficies 
correspondientes a los electrodos se denominan magnodos. 
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La figura muestra un conductor de volumen 
ocupando una región esférica. 

Una bobina enrollada en un núcleo de 
permeabilidad elevada se conecta con el resto 
de un sistema que permite medir la fuerza 
electromotriz inducida entre las terminales a y 
b, e imponer corrientes eléctricas a los 
devanados. Las líneas que salen y entran en los 
extremos del núcleo y atraviesan el conductor 
de volumen corresponden al campo magnético 


que se obtiene al imponer una corriente entre 
las terminales a y b del devanado. 

Las superficies circulares que se encuentran en 
uno y otro extremo del núcleo de la bobina son 
en este caso los magnodos. 


Aplicando la fórmula de reciprocidad al caso que nos interesa, podemos tomar R,=R, 
(siendo R; la región ocupada por la bobina y el resto del sistema de detección) y Ry = Ry 
(siendo Rr la región ocupada por el conductor de volumen de los tejidos) y obtenemos: 


fy, Erplo; PJ (07) aV = f, Eyrlo;7) > J+(w;7) av 


La integral sobre la región R,, de la bobina involucra un producto escalar de dos campos. 

Uno de ellos es la amplitud vectorial de la componente armónica Erj (w; 7) del campo eléctrico 
inducido en la bobina por el campo magnético oscilante producido, en el conductor de volumen 
de los tejidos, por las corrientes eléctricas asociadas a los procesos fisiológicos. 

La amplitud vectorial correspondiente a estas corrientes fisiológicas es] (0; Pda 

El otro campo, J ¿(o; F),es la amplitud vectorial de la componente armónica de la densidad de 
corriente oscilante impuesta a la bobina. 

A través del campo magnético que genera, la corriente impuesta a la bobina induce un campo 


sn 
eléctrico oscilante de amplitud vectorial E), y (w; 7) en el conductor de volumen de los tejidos. 


Dadas las frecuencias involucradas y las dimensiones geométricas de los sistemas considerados, 
la propagación de perturbaciones electromagnéticas se puede despreciar. Se puede trabajar con 
campos casi-estacionarios (Plonsey y Heppner, 1967). 


Plonsey asumió una bobina filiforme. En este caso (Plonsey, 1972) Sa, Es (w; F) Jo; 7) dV 
se puede aproximar por el producto ly(w) V¿(w) . 

En esta aproximación l¿(w) es la amplitud de oscilación armónica de la corriente impuesta a la 
bobina y V,(w) = $ E (0; r(s)) - E (s) ds es la amplitud de oscilación armónica de una fuerza 
electromotriz (fem) inducida en la bobina. 

Esta fem se debe al flujo magnético oscilante que atraviesa los devanados de la bobina y se 
origina en las corrientes eléctricas distribuidas en el conductor de volumen de los tejidos, que 
acompañan la actividad fisiológica del organismo. 
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En la integral de línea $Er» (w; r(s)) -£(s) ds, la función vectorial 7(s) representa la curva 


correspondiente a la bobina filiforme, mientras que t(s) es el vector tangente en el punto de 


1 => Ñ E Ñ 
I(w) Sy Ep,r(w;F) " J7(w; 7) dv 


Reconstruyendo la señal de voltaje en el equipo detector: 
Up (t) = ¡pe expli t w]V,(w) dw 


abscisa curvilínea s. Entonces: V,(w) = 


Plonsey supone: 
(1) Que se puede tomar /¿(w) = 1 para toda frecuencia angular. 


(2) Que Erp(o; Y) es independiente de la frecuencia en el intervalo de frecuencias donde 
Jo; Y) no es nula * y es igual a un campo denominado campo de derivación magnética 
Epu(P). 

En ese caso obtiene: Up(t) = Ja EME) + J ¿(te 7) av 


Se puede observar que para un campo de densidad de corriente JAR) que se mantiene 
constante en el tiempo, en principio la integral podría no ser nula, aunque por supuesto sería 
constante. Pero no debería aparecer una fuerza electromotriz inducida en la bobina cuando las 
corrientes fisiológicas que le dan origen no varían con el curso del tiempo. 

El campo de derivación, independiente siempre del tiempo, se origina en corrientes impuestas 
en la bobina, por lo que no tiene nada que ver con las corrientes que se observan durante la 
actividad fisiológica de los tejidos. 

Cuando la bobina actúa como generadora de un campo magnético oscilante, induce un campo 
de densidades de corriente oscilante en los tejidos, que produce, a su vez, su propio campo 
magnético oscilante. 

Ese campo magnético secundario no se encuentra acotado a la región R7 . 

Su flujo variable puede alcanzar la bobina e inducir allí una corriente que se agregaría a la 
corriente impuesta e influiría a su vez sobre el campo de derivación que aparece en la fórmula 


para up(t). 


Ahora nos proponemos mostrar que: 

(1)-A las frecuencias involucradas en las señales electromagnéticas asociadas a los procesos 
fisiológicos, el campo magnético secundario puede no ser tenido en cuenta cuando el interés se 
focaliza en el cálculo del campo de derivación magnética Enul) producido en los tejidos por la 
corriente impuesta en la bobina. 

(2)-Se puede deducir una fórmula para uv, (t) análoga a la deducida por Plonsey sin necesidad 
de plantear las suposiciones utilizadas por este autor en su deducción. 


20 Lo cual equivale a que la corriente en la bobina es una delta de Dirac. En ese caso la suposición de campos casi estacionarios (o 
sea la posibilidad de no tener en cuenta el tiempo de propagación de las perturbaciones electromagnéticas) es difícil de justificar, 
aún para dimensiones de decenas de centímetros que estamos considerando. 


21 Justifica esta última hipótesis en el comportamiento del SQUID, que en el año (1971) en el que envió su trabajo para ser publicado, 
era un dispositivo de reciente empleo en investigaciones médicas. 
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La corriente de desplazamiento puede ser ignorada para oscilaciones cuya frecuencia angular 


verifica = < 1 donde e y son valores representativos de la permitividad dieléctrica y la 


conductividad del medio. 

Nos proponemos mostrar que si £ es una longitud característica de la parte del conductor de 
volumen donde las corrientes inducidas por la bobina no son despreciables, entonces la 
condición que permite despreciar el campo magnético secundario se reduce a £? wpa <1 


ade rá > a Ñ 
Comencemos con el campo magnético H,,(t, 1) generado cuando se impone una corriente en la 
bobina. 


En el exterior de la región R, y en particular en Ry, se verifica: V X H, (6) 0 


El campo magnético en los tejidos H(t, Y) es la suma de H,(t, F) con un campo secundario 
H.(t, F) debido a las corrientes inducidas. 


Si Er(t, Y) es el campo eléctrico inducido en los tejidos: 
Vx H(t,7)=0 Ep(t,7) Vx Er(t,F) ==10 7; HF) 


Aquí consideramos un tejido homogéneo e isótropo para simplificar. Las conclusiones se aplican 
también al caso general. 

Puesto que difieren muy poco entre sí, se tomó la permeabilidad magnética del vacío como 
permeabilidad de los tejidos en masa. 


=> > 0 7 > 
Entonces: Vx (v xHu(t, 7) == H(t,r) 
Considerando un campo oscilante monocromático e introduciendo las nuevas variables de 


. .s hu y: . 1 . . 
posición Y = £ Y, resulta, siendo V = > Vo y definiendo e=f?4wpu,0: 


Vo Xx (vo Xx É(o, 7.) =-Í€ E(o, To) 


: ey? 
setiene: e= E) 
w Uy O lo) 


Para £=0.25m ?, w = 6.28 x 10* Hz, uy = 1.257 x 102 o y 0<5%/m (una cota 


Introduciendo un espesor de efecto pelicular 9? = 


superior muy elevada para la conductividad de los tejidos en masa) se obtiene Y > 2,5m. 
Entonces € < 0.01 resulta pequeño respecto de 1. 


Desarrollando el campo vectorial fasor en potencias de e€ : 
H(w,To) = HO (w, Fo) + € HV(w, 7) + € HD(w, Fo) + + 


> > 
Sustituyendo en la ecuación Vy X (Vo xH(w, 0) = —i e H(w,F,) se obtiene una secuencia 
de aproximaciones: 
EN 


Vo Xx (vo Xx 40 (o, 7.) =0 Vo Xx (Vo Xx 40D (o, 7.) = 40 (w, ro) 
Vo Xx (vo Xx 4 (o, 7.) == 40D (o, To) ... 


2 £=0.25 m es el radio de una esfera centrada en el corazón, en cuyo interior se encuentran las densidades de corriente más 
intensas generadas por la actividad eléctrica del miocardio. En el caso de la estimulación magnética del cerebro, la dimensión 
característica es menor. 
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De lo dicho previamente sobre el campo directamente generado por la bobina y el campo 
secundario generado por las corrientes inducidas, se desprende que: 


H(w,7o) = Hp(w,Fo) + Hs(w, Fo) 


La amplitud fasorial del campo secundario verifica V¿ X E.(o, Yo) =*2 0 Erlw, F) siendo en 


general Vy X Ér(w,F) = —10 Ho £ É(w,F) 40. 


Esto último junto con V Xx H,(w,7) = 0 nos permite inferir que HO(w,%,) = H,(w,F) es el 
campo vectorial fasor directamente generado por la bobina, mientras que el campo vectorial 
fasor del campo secundario viene dado por: 


H.(w,To) = 4H Vw, ) + 2H (wm, ro) ++ 


Comoe < 1 ,este desarrollo teórico sugiere que podemos despreciar el campo secundario en 
una primera aproximación. Desde el punto de vista práctico un campo secundario tan pequeño 
se confunde con el nivel de ruido magnético. 


Finalmente, reconsideremos la fórmula para la señal de voltaje registrada por la bobina y debida 
al campo magnético generado por los tejidos. 

Es conveniente comenzar con la misma fórmula utilizada por Plonsey (pero sin efectuar sus dos 
hipótesis): 


we) = 72 explico Sa, Eortos7) > T (057) dv) doo 


Cuando a la bobina representada por la función vectorial de la longitud de arco *(s) se le impone 
una corriente ¡,(t) ésta produce en la región Ry un potencial vector: 


> O A t(s) ; 
4 p(t,r) = Le EEE] ds] ip(t) 


La componente del campo eléctrico inducido en los tejidos por la corriente en la bobina viene 
dada por: 


Er) == A = OE $ es) = (00) 2 A (5) 
En esta fórmula £(s) es el vector unitario tangente a la curva que representa a la bobina y la 
norma ||? — F(s)|| expresa la distancia entre un punto F(s) de la bobina y un punto F 
perteneciente a la región ocupada por los tejidos. 
El campo total en los tejidos viene dado por la suma del término debido al potencial vector de 
la bobina y un término debido a un potencial escalar (7: 


Ea > 0 > > 
Ep(t,F) = Voy — a p(t,T) 
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Como V : (v x Ett, 7) =V: (6 . EE, 7) = 0 , de esta última ecuación y de la fórmula para 


el campo eléctrico resulta en el interior del conductor de volumen de los tejidos: 


de d a 
V-(6- Vpr(t,7)) =-V- (7 75% s(67) = [Goo] v: (27 7 Áp(5)) 


Esta es una ecuación para el potencial escalar en los tejidos que puede denominarse ecuación 
de Poisson no homogénea cuasi estacionaria. Cuando la conductividad es un escalar que no 
depende del punto considerado en el conductor de volumen: 


v- (27 0 Áy(5)) =20 0V: ÁL(F)=0 


Entonces el potencial escalar satisface una ecuación de Laplace:  Apry(t,F) =0 


En la frontera del conductor de volumen se anula la componente normal de la densidad de 
corriente, porque el ambiente (aire, ropa) se puede asumir no conductor: 


ñ- (6 -En(t,7)) =0 osea ñ-(5 - Vor(t,7)) =-ñ- (276 AL) LO) 


De la ecuación de Poisson no homogénea cuasi estacionaria y de esta condición en la frontera 


se desprende que el potencial escalar se puede escribir así: prp(t,F) = E lp (0) Py (7) 
El campo escalar Dd (F) es solución de la ecuación de Poisson no homogénea y estacionaria: 
V-(6 - VOp(r)) =-V- (27 5 > Ap (5) 
Junto con la condición en la frontera del conductor de volumen de los tejidos: 
ñ- (6 > VOr(t,7)) =-1- (276 Áy(5)) 
Entonces el campo eléctrico total en los tejidos es proporcional a iy (6) : 


Pese E O) (Vo,(7) +21 A,(5)) = E AO) HENO 


Entonces tanto el potencial vector en los tejidos, dominado por la corriente en la bobina, como 


E A , d. 
el potencial escalar en los tejidos, son proporcionales a qc ib (t) 


Tomando la transformada de Fourier del campo eléctrico Er(t, Y) y denominando lp(w) a la 
transformada de Fourier de ¿,(t): 


>» 1 ns >» > 
Ep(w, 7) = | expl—i w t] [ Ey(t,P)) dt = 1 0w l(w) Ey (P) 
TU —o00 
Sustituyendo este resultado en la fórmula para v,(t) obtenemos: 


O = [7 explito]( fa, Eb 7) +0 Jo; 7) dv) du 
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. . Pi > 0 -i > 
Como fo explitw] iw J5+(w;F) du = Sn J 7(t;7) obtenemos la fórmula siguiente, 


análoga pero no equivalente a la obtenida por Plonsey : 


vy(E) = fp, EN(7) > 5 JH) adv 


t 


A 
En esta nueva expresión, el campo de derivación magnética €, (7) viene dado por: 


Epm(P) = Ep(7) = —VOr(F) — 20 ÁAy(7) (FeR) 


Si el campo de densidad de corrientes es estacionario vu, (t) se anula. 


A de 0 >?¡ > A $ qe 
Si en un punto de los tejidos cn J y(t; Y) resulta ortogonal al campo de derivación magnética, 


no hay aporte a la señal captada por la bobina actuando como receptora, independientemente 


0 >? en 
de la magnitud de 7 J PUT) 


0 Pi > 4 E > E 
Siempre que 7 J 7(t; Y) no sea ortogonal al campo de derivación Epy(F) , en igualdad de 


las demás condiciones, donde la magnitud del campo es mayor, el aporte de las densidades de 


corriente a la señal captada por la bobina va a ser mayor. 


Se puede introducir el campo de derivación Epm(7) como una medida de lo que se denomina 


campo de sensibilidad vectorial del equipo de detección. 


Curva de sensibilidad nula 


El esquema muestra una posición de un 
magnetómetro (unipolar de bobina simple) 
respecto de un conductor de volumen 
representado por una esfera. Las líneas del 
campo de derivación magnética aparecen 
esbozadas como sistemas de curvas planas 
concéntricas respecto de una línea que 
representa los puntos donde la sensibilidad es 
nula (donde el campo de derivación magnética 
se anula). Dado uno de estos planos de curvas, 
que son ortogonales a la línea de sensibilidad 
nula, cuanto más alejado del punto de corte del 
plano con la línea se encuentre un punto del 
conductor de volumen, menor es la magnitud 
|| €, (P)I| de la sensibilidad. 


53 


Esbozo de las líneas del campo de derivación 
magnética y las líneas de igual sensibilidad, 
también para un magnetómetro unipolar de 
bobina simple, pero ubicado en otra posición 
respecto del conductor de volumen 
representado ahora por un el espacio ubicado a 
un lado de un plano. La traza de ese plano con el 
plano de la figura es el eje de abscisas (h = 0) 
donde se ha marcado la distancia r entre el eje 
central de la bobina (eje de ordenadas h que 
pasa por el centro de la bobina y es ortogonal al 
plano de esta última), y un punto en el 
conductor de volumen, cuya distancia a la 
superficie del conductor de volumen se 
representa por h.. El eje r =0 (el eje h ) es la 
línea de sensibilidad cero. 

Las líneas concéntricas a este eje representan las 
líneas del campo de derivación magnética. 

La familia de curvas que se van alejando del eje 
de sensibilidad cero y que en un extremo están 
numeradas (100, 200, 300, 400, 500, 1000, ... 
hasta 107) corresponden a líneas en cuyos 
puntos el campo de derivación posee la misma 
magnitud, es decir ¡puntos donde el 
magnetómetro presenta la misma sensibilidad 
escalar. 


Para la posición de la bobina que muestra la última figura, la sensibilidad escalar lA, Cr, 3) 
debida al potencial vector producido por la bobina se puede estimar por la siguiente fórmula 
4r¿r 
h2+(rp+r)? ” 
elípticas completas de primera y segunda especie: 


donde, por definición k? = mientras K(k) y E(k) son, respectivamente, integrales 


23 


[Ay] =C + /ñ7+ (ro + 11 — k2)K 00 — E(l)) 


En esta fórmula las coordenadas del punto vienen dadas por su distancia al eje r y su distancia 
h a la superficie del conductor de volumen. El radio de la bobina es ry y € es una constante 
numérica. 


La fórmula Ap (7) = (2 $3 ds) se puede aplicar a un magnetómetro formado 


por un sistema de bobinas como el bipolar, como el esbozado en la figura siguiente, y para 


23 La fórmula se obtiene a partir de la expresión para el potencial vector inducido por una corriente ¡¿,(t) en los devanados de una 
bobina circular (Nikolski, 1985; Malmivuo y Plonsey, 1995). 
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magnetómetros multipolares más complejos, siempre que las conexiones entre los devanados 
sean las adecuadas. 


L-- j Li Esquema de un magnetómetro bipolar con 
bobinas ubicadas como polos opuestos. 


En otros casos las bobinas se ubican ortogonales entre sí o del mismo lado y paralelas, 


compartiendo un mismo eje para formar un gradiómetro que permite atenuar gradientes 
espaciales no deseados en el campo magnético. 


En la figura que aparece a continuación se comparan el campo de derivación magnética de un 
par de bobinas opuestas (dibujo en el medio) con el campo de derivación eléctrica de un par de 
electrodos opuestos (dibujo de la izquierda), ubicados en un conductor de volumen esférico. 
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El voltaje que aparece entre los electrodos se expresa como una integral de volumen del 


producto escalar de la densidad local de corriente J5c: Y) por el vector local del campo de 
derivación eléctrica (Malmivuo, 2000; Suárez-Antola 2007 (a) y (b)). 


La superposición de ambos campos, que se puede ver en el dibujo de la derecha, pone en 
evidencia la diferencia en la geometría de uno y otro campo, a lo que se añade la diferencia en 
sus magnitudes. 

En los puntos donde el campo de derivación magnética es nulo, la magnitud del campo de 
derivación eléctrica suele ser máximo respecto de sus valores en los puntos próximos. 
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(13) Consideraciones finales 


-Tanto la evidencia experimental disponible como los trabajos teóricos que plantean 
posibles mecanismos de interacción de los campos magnéticos con los sistemas 
biológicos, sugieren que los efectos no térmicos de la exposición prolongada a campos 
débiles pueden tener consecuencias útiles y consecuencias nocivas desde el punto 
biológico y médico. La posibilidad de la existencia de ventanas de baja frecuencia dentro 
de las cuales las interacciones se pueden amplificar por resonancias no puede ser 


descartada. 


-Las normas actualmente en vigencia para limitar las exposiciones a los campos 
eléctricos, magnéticos y electromagnéticos se basan en la investigación de umbrales 
más allá de los cuales se producen efectos agudos (que aparecen a lo sumo en minutos) 
añadiendo luego factores de seguridad conservadores. No obstante, hay evidencia 
acerca de efectos de los campos magnéticos sobre el crecimiento celular, que se 
produce fuera de esa ventana temporal. Por otro lado, procesos más lentos dan tiempo 
para que los mecanismos de reparación de alteraciones estructurales, con 
consecuencias funcionales nocivas, puedan actuar y si lo hacen, ciertos efectos podrían 


no aparecer. 


-En el caso eléctrico lo que se combina con el correspondiente campo de derivación es 
la densidad local de corriente, mientras que en el caso magnético lo que se combina 
con el campo de derivación es la derivada parcial de esa misma densidad de corriente 
respecto del tiempo. A partir del ejemplo de las líneas de los campos de derivación 
eléctrica y de los campos de derivación magnética esbozados en la Figura 16 se 
comprende que, en general, la información que se obtiene midiendo fuerzas 
electromotrices producidas por campos magnéticos variables es complementaria 
respecto de la información que se obtiene midiendo el voltaje entre electrodos para una 
misma distribución de corriente debida a procesos fisiológicos. 
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